NASA CR 120985 
BOEING D38677 


TOPICAL REPORT 


ACOUSTIC ATTENUATION ANALYSIS PROGRAM 
FOR DUCTS WITH MEAN FLOW 


by 

R. K. Kunze, Jr. 

(NASA-CB-1 20985) ACOUSTIC ATTENUATION 
ANALYSIS PROGRAM FOR DUCTS WITH MEAN 
FLOW fftoeing Co., Wichita, Kans.) 95 p 
; * \ CSC! 

' — THEBOEING COMPANY 

3801 SOUTH OLIVER 
WICHITA, KANSAS 67210 


Prepared For 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

August 17, 1972 


CONTRACT NASA 3*14321 


NASA-LEW1S RESEARCH CENTER 
CLEVELAND. OHIO 
H. BLOOMER, PROJECT MANAGER 
V/STOL AND NOISE DIVISION 


( Reproduced by 

' NATIONAL TECHNICAL 

i INFORMATION SERVICE 

1 US Department of Commerce 

J Springfield, VA. 22151 




N75- 11211 

20D Unci as 

G3/3 4 .. 02767 v 



FOREWORD 


Ths work described herein was dona by The Booing Company, 
Wichita Division, under NASA contract NAS 3-14321 with Mr. 
H. Bloomer, V/STOL end Noise Division, NASA • Lewis 
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ABSTRACT 


A computerized acoustic attenuation prediction procedure 
has been developed to evaluate acoustically lined ducts for 
various geometric and environmental parameters. The 
analysis procedure is based on solutions to the acoustic 
wave equation, assuming uniform airflow on a duct cross 
section, combined with appropriate mathematical lining 
impedance models. The impedance models included in the 
analysis procedure are representative of either perforated 
sheet or porous polyimide impregnated fiberglass facing 
sheet coupled with a cellular backing space. Advantages 
and limitations of the analysis procedure are reviewed. 
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NOMENCLATURE 


v- ' 


C 

d 

f 


h X( h 


y 


k 


k x , ky, k z 
k r , m, k z 
n 


P 

t 

r, 0, z 


x, Y, z 


P, q, r, s 

J V , Yy 

H* W.H V (2) 
L 

u 

M 

R 

U 

V 

z 

* ' 
v 

P 

CJ 

M 


Velocity of sound, cm/sec 
Acoustic lining backing depth, cm 
Frequency, Hz 

Rectangular duct cross section dimension, cm 

k = ufc 1/cm 

Rectangular duct wave numbers 

Annular duct wave numbers 

Hard wall mode number 

2 

Acoustic pressure, dynes/cm* 

Time, sec 

Cylindrical coordinates 

Rectangular coordinates 

Cross products of Bessel functions 

Bessel functions, first and second kind, respectively, of order v 
Hankel functions of order v , Type 1 and Type 2 
Lining treatment length, cm 
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Attenuation rate, db/cm, db/inch 
Duct mean airflow velocity, cm/sec 
Acoustic particle velocity, cm/sec 
Specific acoustic impedance c.g.s. rayls 
Particle displacement, cm 
Bessel and Hankel function order 

q 

Density, gm/cnrr 

Angular frequency, radians/sec 

Viscosity (poise) 
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SUMMARY NASA Fan Noise Suppression 

A systematic acoustic lining design procedure was de^ope^^^ jncluding the Acoustic 
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pressure, Mach number, and the acoustic spe ^orated 

The impedance models included in the 
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dependence on material porosity characteristics, temperat 

The basic analysis determines the modal P r °P a f*°" 
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performance of the total ^ch Variations m ^^"'procedureles neglect the effects of 

r -rr: see « »- 
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presented in this report. 
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1.0 


INTRODUCTION 


1.1 Background 

A systems acoustic lining design procedure was developed for the NASA Fan Noise 
Suppression Program. This design procedure consists of several computerized modules, 
including the Acoustic Attenuation Analysis Program. The program combines the basics 
of two technologies, acoustic propagated wave analysis and material acoustic impedance 
modeling, into an acoustic lined duct performance simulation. 

The zero airflow work of Cremer (Reference 1) and the mean flow results of Eversman 
(Reference 2) were incorporated into an early version of this program, using a 
rudimentary impedance model for perforated sheet impedance (Reference 3). This 
version was utilized in the NASA Langley Treated Tailpipe Program (Reference 4). More 
recent studies at Boeing-Wichita (Reference 5) have led to improvements in the 
impedance model for perforated sheet and development of an impedance model for 
polyimide liners. A portion of the experimental impedance data for these models was 
included in work reported by Rice (Reference 6). These models have been incorporated 
along with expanded wave analysis applications, into the current version of the program. 
A detailed description of the analytical techniques is presented in the following sections. 


1.2 Technical Approach 

The technical approach used in the development of this analytical program is the 
solution of the mean flow acoustic duct wave equation. The duct is assumed to be 
infinite in length with one or two walls lined. A rectangular duct is assumed for 
computational purposes, and annular geometry applications are treated by 
approximating a sector of the annulus by a rectangle. Normally reacting duct walls are 
. assumed for lined surfaces. Variations in Mach number, wall impedance, and duct 
dimensions in the direction of airflow are treated approximately, and acoustic wave 
reflections due to these variations are neglected. Wall impedances are provided by 
semiempirical impedance models for perforated plate and polyimide facing sheets with 
cellular air backing cavities. 

The computer application of the acoustic wave equation solution to the- acoustic lined 
duct performance simulation is shown in the block diagram of Figure (1). The major 
' components of this program will be discussed in detail in the following sections. 
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2.0 


ANALYTICAL MODEL 


The analytical model for this study is that of a semi-infinite duct with airflow. The 
airflow is assumed inviscid and nonturbulent, with a uniform axial velocity U on a duct 
cross section. The acoustic pressure is assumed to be small in magnitude in comparison 
with the uniform pressure. 


For a duct with the flow axis in the Z coordinate direction, the acoustic wave equation 
is 


C 2 v 2 P = 



( 1 ) 


where 


c = velocity of sound 


p - acoustic pressure 
U = uniform airflow velocity 


The solution and boundary condition analysis for a rectangular duct, as depicted in 
Figure 2, follows. 


An assumed solution to ( 1 ) is 


P = e icot e~' ik z 2 j^Ae ik x x + B e ~ ik x*j Ce'M + De- |k yyj 


( 2 ) 


Substitution of this solution into (1) results in the following relationship between the 
wave numbers k, k x , ky and k z ; 

k 2 = (_-L^)(-kMi 1 /k2 -(l-M 2 )(k x 2 . k y 2)) • -(3) 


where 


to - angular frequency - radians/sec 
k = o>/c 
M - U/c 

The expressions for k x and ky are determined from boundary conditions at the duct 
walls. These conditions are obtained from the assumption of continuity of particle 
displacement at the walls. Let ^ be the particle displacement normal to the wall at X = 
h x . The equation of motion for fluid flow at X = h x in this direction is 



where 
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p = air density gm/cm^ and (D/Dt) is the total derivative with 
respect to time (t) 

D _ _9_ 9_x P_y _9_ 9z_ _9_ 

Dt " Pt + Pt 5x + Py Pt Pz 

For the wall at x - 0, (4) becomes 

^(£<M 2 r (SrL 0 !51 

where £2 * s the particle displacement normal to the wall at x « 0. The difference in sign 
between (4) and (5) is due to the change in direction of the normal to the wall with 
respect to the x coordinate axis. 

The mean velocity components in the x and y direction are zero, thus (4) and (5) 
become, respectively 



The acoustic impedance (Z/pc) and admittance { L ) may be defined by 


I P_ 

pt 


Z 

PC 



(8) 


I 

PC 


P 

££ 

Pt 



(91 


For harmonic motion 


P 

Pt 


tj) = i 


a 


% 


and from (2) 


P 

91 


P 


- ik< 


p 
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applying these relations to (6), (7), (8) and (9) 


we obtain the boundary value equations 



<l-Mk z /k) 2 P 


(l-Mk 2 /k) 2 p 



PP 


x = o 


( 10 ) 


( 11 ) 


application of the boundary condition (10) and (11) 
h x leads to 


to a duct with one lined waH at X = 


at y 2 01 1 = 0 and = o 

£ y 


y = h y : L = 0 and — = o 

<?y 


x = 0: L = 0 and = o 


at x - h x : L = L, 


and 



for the y direction we obtain 
SIN (kyhy) = o 



0 , 1 , 2 , . 


( 12 ) 
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for the x direction 


iLjkfl COS (kxhx) = -kxSIN(k x h x ) 

or X k ' 

iLjkhx = hxkx TAN {h x k x )/^l - 

In the case of a duct with equal admittance walls at both x = 0 and x = h x , we obtaii 
x boundary equations in the following manner from (2), ( 1 0), and (11)/ 

, L|k ( l .“Ji) 2 [ A ,»«*hx * Be-*xh*] 

= -ikxfAe'Mx . Be' ik x h *J 


iL '"( l -nr) 2 [ A * B ] 

- ikx j^A - bJ 


where L-j is the admittance at x = 0 and x - h x 

These two relations may be combined to eliminate A and B, resulting in a qua 

(iL,kh x ): 


h x 2 (iL|k) 2 (l- -^) 4 SIN(h x k x ) 

/ Mk z \ 2 o o 

+ 2iLjh x k ^ I - ■ ~ — J h x k x COS k x h x - h x ^^ x S IN k x h x 


= 0 


or 

-hxk x (COS(hxk x )tl)/SIN(kxhx) 
(I^Mk z /k) 2 

«6 


iL|h x k = 



Wng half angle trigonometric 


for the 


identities (17) becomes 


Sign 


tL|kft x / 2 = _.^*kx /2)C0T (h x |< x /gj 

(I^Mk7A)2 


for thp « 
me - Slgn 


iL|kh x /2 - (hxk>/p) 

(< - Mk z /k)2 



3 -0 METHOD OF APPLICATION 

The acoustic attenuation prediction program provides an analytical evaluation of the 
attenuation characteristics of an acoustically treated rectangular duct. Lining 
configurations may be either a single treated wall or two opposite walls with identical 

r ~' n add 1 ,t,on ' P r °P a 9at*on of the acoustic pressure wave may be in opposition 
to the flow or in the same direction as the flow. 

3.1 General Description 

The analytical acoustic attenuation prediction program evaluates the attenuation 
c aracteristics of rectangular ducts with mean airflow and acoustically lined walls, having 
either one or two walls treated. a 

The basic analysis consists of determining the acoustic modal propagation constants for a 
given duct geometry, environment, and lined wall impedance. The environmental 
parameters considered are temperature, pressure, Mach number, and duct acoustic source 

spectrum The lined wall impedance is obtained from acoustic impedance mathematical 
models, which is discussed in Section 3.2 

Initially, the input acoustic source spectrum is subdivided at each frequency into modal 
pressures. The wall impedance at each frequency is then determined for the first lining 
increment, and the boundary value problem solutions are obtained. The resulting modal 

attenuations are now applied to the modal pressures, and the corresponding attenuated 
acoustic spectrum is calculated. 

' S rePea ! ed for each linin 9 '"cement until the entire lined length has been 
analyzed. However, the source spectrum is subdivided into modal pressures for the first 
increment only, with each succeeding incremental treatment utilizing the attenuated 
modal pressure spectrum resulting from the preceding increment. Also, the resultant 

acoustic spectrum at each increment is used for impedance determination in the next 
increment. 


3.2 impedance Models 

The imp e dan ce mo de|s in the present program include mathematical representations for 
-both perforated plate and porous polyimide facing sheets with cellular air backings 
terminated by an impervious sheet {Figure 3). ® 

I h - m0de ' S COnSiSt ° f a mathematical description of the impedance of the porous face 
sheet dependence on material porosity characteristics, temperature, pressure and total 
particle velocity The total particle velocity V T is a combination 'of S vefoch e 

' e e xchation V m ^ 9raZi " 9 V the face a " d -ustfc pr^um 


a - 

V T = ( Vgf 2 + v a 2) f? BCa)Ii,G PAG£ NOT FILMED! 


( 20 ) 
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TK. „p« «P»-« <« 

perforated sheet, including expressions for the respective gra y 
are as follows: 


Polyimide 
Z = Q + .305 P N 


2, + 2.86 PN 2 Vp + i { .549 P N 3/2 f • 5.12 -lO' 6 P N 3 1Vp} (21) 


.05cM 2 

V ,= — — i ’ 

9 [q/Pc + •J (Q/ pc) 2 + .71 N 2 M 2 

where Q - 4640. P N 

i = viscosity (Poise) 

i =yr 

c - velocity of sound in air (cm/SEC) 


N - number of plies of polyimide 
P = density of air (gm/cm ) 

Vp = particle velocity 
f = frequency 
Perforated Sheet 


( 22 ) 


Z/p c * Q + 


3.83-1O- 5 0 -75VT 

|t/D + 1 - P 0A 1 + tt/2 VT Vp/c Ee * 77 

P 0A yfdleVM^T L J 


4.69- 10~ 4 f 


(T + d)} 


.25 KcM' 


= 


gf Q/2 +n/(Q/2) 2 + K/2 EM 2 
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where 


These 




r^;ioj v P 




= .85D (1 - -7v 


f = frequency (H z ) 

2.54 D f PoA 
St > = ^ 

T = sheet thickness (inches) (cm.) 

hole diameter (inches) (cm.) 

o D wni9(=T T (°K)/288.33) 

* ■ T r"“ u ::; P »T 

, iHii =f — 1 


q = perforate 


<5 = 


,077 T / _£1 \ 

° = ~^? oa \6 + -416/ 

E = 1.0251 (VP 0 A 2 - 1H1/P ° A> 
K = .05 + .11 D/<9rj 

erj = boundary layer momentum 

\/p - particle velocity 


■1 e - .5072 T/D 


thickness (inches) (cm.) 


face sheet impedance 


models are described in greater 


detail in Reference (5). 
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wall impedance, Z w , is given by 

7 » Zc - i cot (kd) 

*-w r 

where d is the depth of the air backing cavity. 

The acoustic pressure particle velocity due to the acoustic pressure at a smgle frequen y 
Pj is given by ^ 

v ai " Z w 

Particle velocities and ^ 


In the first case, the expression for total particle velocity calculation is 


V T {f) = (v gf 2 


1/2 


(26) 


where 


f = frequency 

V T {f) = total particle velocity at frequency f 
P(f) = acoustic pressure at frequency f 

z (f v T f) - wall impedance at frequency f and particle velocity V T (f) 

**w. ' • , 

In the second case, the total particle velocity is denoted to indicate 

• j xUn Avm*Dccinn hpcomes 


(27) 


In tne sutuiiu — 

dependence on acoustic pressure at all frequencies, and the expression becomes 

v T „, - (V *2 V) 

Pf 

where V o f = ]Zw f < v T rms '| 

. n , are solved by iteration techniques based on the Newton Raphson 
. These expressions are solved oy impedances. 

method to obtain V T (f> or V T and the wan 
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3.3 


Modal Pressure Subdivision and Recombination 


The acoustic pressure at a given frequency is subdivided according to the number of 
modes which can propagate at that frequency. The cutoff frequency for any given mode 

number n is given by f = J l g f° r unlined ducts. The cutoff frequency is 

decreased for lined ducts, and a factor of 1.1 has been arbitrarily introduced in the 
denominator of the above equation to account for this effect. 

The subdivision of acoustic pressure into modal pressure amplitudes is accomplished 
assuming the acoustic modes are uncorrelated and initially have equal modal pressure 
amplitudes. The total sound pressure level at any frequency ($PL tota |) is given in terms 
of the total acoustic pressure (P^otal^ ^y 


SPLiolal * 10 logi 0 (Ptota|/ 00022) 2 (28) 

where Ptotal ' s ' n dynes/cm 2 

Defining SPL: as the sound pressure level of the ith modal pressure amplitude, for 

'modal 

n propagating modes the expression for modal amplitude subdivision is 


SPL 'modal = SPL »otal - I°9|0 n 


(29) 


The recombination of attenuated modal sound pressure levels to obtain a total 
attenuated SPL is given by 


SPL total = 10 ,0 9 10 


2 IO ( SPLi modal/lo) 


(30) 


where all SPL's in this expression have been attenuated. The attenuation for any 
frequency is then the difference between the initial and attenuated values of SPL totai . 


3.4 Boundary Value Solution Application 

The boundary value problem solutions are obtained for each frequency of interest and 
for each increment of wall lining analyzed. These solutions are obtained by iterative 
procedures based on the Newton Raphson method. Rewriting the boundary value 
equation for the symmetric modes of a duct with two walls lined, equation (19) becomes 

iL j kh x / 2 = (h x /2)k x tan (h x /2)k x /(I - Mk z /k) 2 (31) 
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It Is noted that the required solutions for physical wall impedances lie in the first 
quadrant of the complex plane {h x /2)k x . {The third quadrant is a mirror image of the 
first, while the second and fourth quadrant corresponds to wall impedances whose real 
part is negative.) 

The solution to (31) for zero Mach number and unlined walls are easily obtained real 
quantities. 

The starting points for the iterative solutions for no-flow, lined ducts are chosen near the 
unlined wall solutions. Similarly, the lined wall duct with flow is solved using the lined 
wall no-flow solutions as starting points. However, the interative solution for the flow 
case must be accomplished in incremental Mach number steps to obtain all solutions of 
interest. The incrementation of Mach number by .05 steps appears to be sufficient for 
most applications. 

In some cases, the iterative procedure crosses into the second or fourth quadrant of the 
complex plane. This effect is checked in the solution procedure, and countered by 
returning the iteration to a new starting point. This point is specified by selecting the 
point on a line segment connecting the non-first quadrant point to the last iteration point 
in the first quadrant. 


Modal Ordering 

The eigenvalue solutions are ordered according to increasing modal attenuation rate. This 
ordering scheme assures that the least attenuated mode is always considered. Under the 
above stated assumption of initially equal modal pressure amplitude, it can be seen that 
the least attenuated mode is the most significant contributor to the resultant sound 
pressure level after some length L of treatment. Consider the simplified case of two 
modes with attenuation rates R-j and R 2 . If the final total sound pressure level is 


SPLt = 10 log, 0 (l0SPL|/l0 + , 0 SPL2/10) 


(32) 


where SPL-j and SP1_ 2 are the sound pressure levels of modes 1 and 2, respectively, after 
length L. The difference, A , between SPLy and SPL^ can be expressed as: 


A = 10 log (0 (I + I0-< r 2 “ R |) L/, °) 


{33) 


From this representation, it can be seen that the maximum difference is 3 dB, which 
results if R j - R 2 and that for R 2 “ R i + *58, the difference becomes 1 dB. 


Lining Length Incrementation 

The analysis of a duct by incremental lengths is incorporated to account for effects of 



changes in Mach number, duct, height, and lining characteristics with progressive 
position in the direction of acoustic propagation. These effects are treated by 
specification of the varying parameters for each increment analyzed. Acoustic wave 
reflection is not considered, and thus the analysis is applicable for small parameter 
variations. 
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4,0 PROGRAM APPLICATION AND LIMITATION 

The assumptions of acoustic modal content, mean airflow, and neglection of acoustic 
wave reflections are the major limiting factors in the application of this program. Some 
of the pertinent characteristics of these assumptions as related to program application 
are as follows: 


4.1 Modal Content 

The equal modal pressure amplitude assumption is incorporated since experimental 
modal pressure amplitude definition is unavailable, and for ease of application. As 
definitive model content information becomes available, it can be incorporated in this 
program. 

4.2 Duct Length Incrementation 

The sensitivity of the impedance models to sound pressure level spectrum variation 
requires the analysis of an acoustic liner in several incremental lengths. A check on the 
validity of any increment length utilized may be accomplished by comparison of the 
attenuation rates of each lining segment; gross differences indicate the need for further 
subdivision of the increment length used. 

A further consideration for increment selection is the degree of variation of the other 
environmental parameters with length, e.g., Mach number, geometry, etc. Variations in 
these parameters may require a further reduction in the increment length than for fixed 
parameters. 

4.3 Multiple Lining Analysis 

The treatment of cases having sequential segments of differing acoustic liners should be 
accomplished with one program execution in preference to an execution for each lining 
type. This approach maintains the continuity of modal pressure contributions for the 
entire lining treatment, conversely, an execution for each lining type would redistribute 
^ the modal pressures, and the final duct attenuations would be overestimated. 

4.4 Application to Annular Geometry Ducts 

The subject program can be utilized in the performance prediction of annular geometry 
ducts by approximating a section of the duct annulus with a rectangle. Details of this 
, type of application are treated in Appendix I. 

4.5 Application to Ducts with a Sheared Flow 

The applicability of the subject program to ducts with a sheared flow (nonuniform 
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velocity profile on a duct cross section) depends upon the relative direction of the 
acoustic wave propagation and the airflow, and upon the magnitude of the boundary 
layer thickness and Mach number. Using the nomenclature “inlet mode" and "exhaust 
mode" to denote cases in which the acoustic flow and airflow are In the opposing and 
identical direction, respectively, the work of Eversman (Reference 7) and Munges and 
Plumblee (Reference 8) indicate that the shear flow effects are generally small for the 
exhaust mode, but can be very significant in the inlet mode. Analytical methods have 
been developed to correct for the sheared flow effect, but have not been included in this 
document. 
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5.0 


CONCLUDING REMARKS 


The acoustic attenuation analysis program provides an analytical capability to evaluate 
lined duct configurations having uniform flow on a duct cross section. Variations in duct 
geometry, Mach number, and lining impedance are permissible within the duct. However, 
variations in these parameters are treated approximately in that acoustic wave reflection 
effects are neglected. Initial modal pressure amplitudes are assumed equal, modal 
interaction effects are neglected and modes are assumed to be in phase. Acoustic wall 
impedances are obtained from semiempirical impedance models. 


The acoustic attenuation analysis program was developed primarily for application to 
aircraft fan jet engines. The major limitation is considered to be the assumption of mean 
airflow within the inlets of these engines. For this case, the program results should be 

corrected to account for the sheared flow effects to obtain a more accurate estimate of 
the acoustic attenuation. 

Further analysis and subsequent program improvements are indicated in the followinq 
areas: 


• Sheared flow effects 

• Modal pressure amplitude definition 

• Modal phasing and interaction 

• Acoustic wave reflection effects 

The first item requires a reformulation of the analytical procedure so that the sheared 
flow effects can be evaluated as an integral part of the analysis program. This 
modification would provide a more detailed and rigorous analysis than correcting 
the mean flow results. 


The second and third items are contingent on having adequate space-time data from which 
modal amplitude and phase information can be extracted. These data are necessary to 
compute the initial modal energies and to track these energies down the treated duct. 

The fourth item deals with the reflection of acoustic modes in the treated duct. These 
reflections can result from changes in the wall impedance and/or geometry. This problem 

can be approached using either assumed modal amplitudes and phase or when available 
the results of items two and three. 
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INPUT 

PARAMETERS 


INITIAL MODAL 
PRESSURE AMPLITUDE 
DETERMINATION 


LINING INCREMENTATION 
/UPDATE MACH NO,\ 

(geometry, lining) 

V SPECIFICATIONS / 


FREQUENCY 

INCREMENTATION 


LINING 

IMPEDANCE 


WAVE EQUATION 
BOUNDARY VALUE 
SOLUTIONS 
AND MODAL 
ATTENUATIONS 


ORDER MODAL 
SOLUTIONS 


MODAL ATTENUATION 
RECOMBINATION 
AND TOTAL ATTENUATION 
DETERMINATION 


OUTPUT 

ANALYSIS PROCEDURE FLOWCHART 
FIGURE 1 
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RECTANGULAR DUCT GEOMETRY 
FIGURE 2 
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BACKING 
DEPTH (d) 



ACOUSTIC LINING PANEL 
FIGURE 3 
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format 


VARIABLE NAMES 


72 Column Alphanumeric Field 

8 Fields, 8 Columns each 

72 Column Alphanumeric Field 

9 Fields, 8 Columns each 

« if 


Title 

A, SWM, FWM, H, PSPL, RMS, REORD, OVRMOD 
SPL NAM 

SPL 0), SPL (2) SPL (9) 

SPL {10), - 


SPL {99} 

TL, DL, DHH,G,TT, PT, A, B 

ANPLY, AM, DHS, OPA NP, D, DIA, R, MTHK 

AAM (1), A AM {2) 


ADHS (1), ADHS (2) 


P0AC{1), POAC {2) 


DCON {1), DCON {2) . 


DIAC {1), DIAC (2) 


THCK {1), THCK {2) 


MTHK {1), MTHK (2) 


INPUT FORMAT 
FIGURE 5 
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APPENDIX I 


ANNULAR GEOMETRY ANALYTICAL MODEL, APPROXIMATION 
BY RECTANGULAR GEOMETRY 


The annular geometry analysis follows that of the rectangular duct solution, using cylindrical 
coordinates in the wave equation solution. The duct geometry and coordinate system is shown in 
Figure 4. The acoustic wave equation is again 

c2v2p = (rf + u §if p <A1) 

Separation of variables leads to the solution: 


P = (AJ m (rk r ) + BY m {rk r ))cos(m'$)ei( Cu t - k 2 Z) 
and 


(A-2) 


k z = (irW) (' Mk " /k 2 -(l-M 2 )k r 2 ^ 


(A-3) 


where r, 0 , z are the system cylindrical coordinates 
k = cj/c 
M = U/c 

m - angular wave number = 0, 1, 2, 

■W rk r> - Bessel function, first kind, order m 
Y m< rk r> “ Bessell function, second kind, order m 

Both walls are assumed lined with identical wall admittance I". The boundary condition at r , the 
outer annular wall, and rj, the inner annular wail for this assumption become 


at r i • iKL(,-M-^) 2 P 
at r 0 : ikL (l - M ^j 2 P = - ^ 


Substitution of A*2 into A-4 and A-5 and solving, the following eigenvalue equation is obtained: 

PRWWING PAGE BLANK NOT 
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- 2 g 

j"(iKL)(l - M-j -) 2 p m ' t'kL) (> - M ~ ") (<lm - r m> ~ S m r 0 


(A-6) 


where 


a = r 0 k r 


b = rjk r 


fy - J>»Yy(b) - Jy(b)Yy (a) 

q v = (jytaJY^'lb) - jy(b)Yy(a))k r 
ry = (j/laJY^b) - J/b) y/ ( a)) k r 
= (j/(o)Y/(b) - j/(b)Yy(a)}k r 2 


(A-7) 

(A-8) 

(A-9) 

(A-10) 


(Primes denote differentiation of the Bessel function with respect to the particular function 
argument.) 

As in the rectangular case, A-6 may be written as the quadratic solution 


ikL z ~2p ” r m ^ ~ r nrT 




,a "> 


The approximation of an annulus by a rectangle might be presumed, by inspection, to be valid 
whenever the radius ratio, r-./r., is close to unity. This presumption may be mathematically verified 
in part by applying large argument Bessel function approximations to A-6. Taking the definition for 
Hankel functions in terms of Bessel functions: 

H-y (l) (Z) = j y(Z) ♦ iY/(Z) 


H / (2)(Z) = J r (Z) - i Y r (Z) 

where Hy W(Z) and H y * 2 *(Z) are Hankel functions, type 1 and type 2, respectively. The large 
argument approximations for H y ^(ZlandHy (2) (Z) are (Reference 9): 


H^t'hz) ~ JZ/TrZ e'^-^-'-i 1 ) 


(A-12) 
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(A-13) 


Hy (2) (Z)~ fz/7rz e -i(z--^£--^.) 

H a n kef fun cti o n c^oss product ' W ^ ^ ' A " ? ' A ' 8 ' A ' 9 and A " 10 ™ V be Wri,ten in terms of 

p y ' - ■jr(H^< | )(a)H/-<2)(b) - < 1 ) (b) H r < 2 )(a)) (A-14) 

q y ' ( 1 ^(a) (2) (b) - O’tblHy ( 2 )(a» (A-15) 

V = ■'57 (H / (I) (°) (2 >(t>) - H ^«)( b )H r (2)'( 0 )) (A . 16) 

k 2 i , 

S y 1 * ~ir (H r (l) (o>Hy< 2 >'(b) - (l)'(b) H/ (2)'(o)) (A-1 7) 

Applying A- 12 and A-13 to the above yields: 


*Y~ 7TJ Hb S ' N(a - b) 


(A-1 8) 

q > ~ k '(^ cos(a - b) + T p x) 


(A-1 9) 

r / ~ k 'Ua b cos <°-» -f p y) 


(A-20) 

^ " - kf2 ((ab +l )( W ) SIN! °- b) + 1 

(f-f)^b C0S( °- b) ) 

(A-21) 


Using the above relations in A-1 1 results in 

' L|k_ 2SIN(o-b) {- 2COS < a - b ) -^)siN(o-b) i (A-22) 

b') S,N2(a - b)t8 (7--^) c 0S(a-b)SIN{a-b)]^2j 


ioting (a - b) h k r where h is the distance between lined walls for the 
/vhere r = o (radial modes only), A-22 becomes 


annular duct. 


For the case 
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( A-23) 


i L jkh = hkr 


(COS(hkr) t I) //, 
SIN(hkr) / \ 



which is identically the rectangular solution (Equation 17). 

The use of the large argument approximation has two implications', first, for large radii and/or kj, 
the rectangular approximation appears valid, and second, a further limitation on the relative 
magnitudes of r 0 and Tj is implied. This limitation is that if r Q k r satisfies the large argument 
approximation, then r i /r Q should not be nearly zero, since the approximation is applied to both 
r Q k r and rj k r 
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APPENDIX II 


program description 

This program consists of a mainline control program and subroutines. The mainline program calls 
subroutine 'input" to obtain all input data and obtain preliminary processed data. The 
incrementation loop is then set up and the subroutine "ATNSPL" is called for each lining 
increment. "ATNSPL" serves as overall control for ail calculations for each increment except for 
the calculation of an attenuation spectra at the current increment. A maximum of 50 increments 
may be treated. 


Subroutine "INPUT" 

This subroutine reads all input data, calculates speed of sound, density, and the numbers of 
propagating modes for all frequencies, and model pressure spectra based on these data. The input 
format is shown in Figure 5. The variable name definitions and the input options are as follows: 

Identification only 

Number of 100 Hz bandwidth SPL values to be input (beginning with 
frequency of 200 Hz, maximum value 99) 

Maximum number of soft wall modes treated 

Maximum number of hard wad modes treated 

Number of walls lined (1 of 2) 

Flat SPL spectrum value; if zero, will read "A" values of SPL, if nonzero, 
sets up flat spectrum array of "A" values of SPL equal to "PSPL" 

0. Impedance programs calculate particle velocity and impedance as 
function of SPL at one frequency 

1. Impedance programs calculate an "RMS" particle velocity based on 
entire SPL spectrum 

REORD 0. Modal ordering by least attenuated mode 

1. Modal ordering by largest real part of downstream wave numbers 
OV RMOD 0. T reats number of modes propagating to limits of "SWM", "HWM" 

1. Treats "SWM", "HWM" modes, cutoff conditions not applied 
SPL NAM SPL spectrum identification title 


TITLE 

A 

SWM 

HWM 

H 

PSPL 

RMS 
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SPL Array of up to “A" sound pressure levels 

Tl Total duct lined length (inches) 

DL Increment length (Inches) 

DHH Hard wall duct height (inches) 

G Ratio of specific heats for duct environment 

TT Total temperature (°R) 

PT Total pressure (PSIA) 

"A", "B" Frequency limits for calculation of attenuation (H z ) 
ANPLY Impedance model selection 
0. Perforated sheet 


1, Polyimide 
AM Mach number 

DHS Lined wall duct height (inches) 

OPA NP Perforated sheet fractional open area or polyimide number of plies 

D Backing depth (inches) 

DIA Hole diameter (inches) 

R Material thickness (inches) 

MTHK Boundary layer momentum thickness (inches) 

("DIA", "R", and "MTHK" are required for perforated sheet impedance model 
only.) 

The following data arrays are required when the corresponding variables above are zero. The 
number of values per array is equal to the number of increments as determined from the total 
length and incremental length: 

AAM(I) Mach number 

ADHS(I) Lined wall duct height (inches) 
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POAC(I) 


Perforated sheet fractional open area or number of polyimide plies 
DCQN(l) Liner backing depth (inches) 

DIAS(I) Perforated sheet hole diameter (inches) 

THKC(l) Perforated sheet thickness (inches) 

MTHKC(I) Boundary layer momentum thickness (inches) 

Subroutine "ATNSPL" 


Subroutine "ATNSPL" is the control routine for the calculations pertaining to one duct increment. 
Initially, if RMS impedance treatment has been selected, the RMS particle velocity for the desired 
lining type is obtained by a call sequence to the "ENTRY" section of the appropriate impedance 
subroutine. If the discrete frequency impedance treatment option has been selected, the above call 
sequence is skipped over, and the frequency increment "do loop” is entered. 

The frequency increment "do loop" controls calculation of impedance, eigenvalue equation 
solutions, modal attenuation ordering, and attenuated SPL spectra for each frequency in the range 
specified. 

The impedance is obtained by standard calls to one of the impedance routines "IMPQD" or 
"IMPPD". A parameter "TSP" in the impedance routine call sequence specifies whether the routine 
is to calculate discrete or RMS particle velocity impedance. The impedance and other required data 
are then entered to the subroutine "UX" which determines eigenvalue solutions for the lined walls. 

This subroutine returns the solutions ordered in terms of increasing attenuation. The subroutine 
"RORDER" is then called if ordering is desired according to decreasing real part of the down 
stream propagation constant. 

The (n,o) modal attenuations are next applied to the appropriate modal spectra. If more than one 
hard wall mode is specified, the remaining eigenvalue solutions are obtained from subroutine 
"HWDBS" and the resultant attenuations applied to the modal spectra. Finally, the attenuated SPL 
spectra are synthesized from the modal spectra, and control returned to the "Mainline" program. 


Impedance Subroutines "IMPQD" and "IMPPD" 


These subroutines are the impedance models of the program. The models are semiempirical models 
for perforated sheet and polyimide, respectively, with cellular air backing. 

Each routine is set up in two parts — the first accessed by a standard subroutine call, and the second 
accessed by an "entry" call. The second part is used to determine an "RMS" particle velocity for 
the lining, while the first part calculates either an impedance from an imput "RMS" particle 
velocity or an impedance based on a discrete frequency particle velocity. 
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Subroutine "UX" 


This subroutine controls the calculation of lined wall eigenvalue solutions for one frequency. Th; 
zero Mach number solutions are obtained from subroutine "RTZ" and are filtered for duplicate am 
invalid solutions. The filtered set of solutions are then input to subroutine "MRT" from whicf 
solutions for the required Mach number are obtained. Control is then returned to subroutim 
"ATNSPL". 


Subroutines "RTZ" and "NRAPH" 

"RTZ" controls the calculation of zero Mach number solutions accomplished in "NRAPH". Thr 
maximum number of modes possible is limited only by the size of the array "ALZRT" in "RTZ' 
and "UX", and array "A" in "ANTSPL". Presently, these arrays are dimensioned at 42, which give: 
a maximum of twenty modes. If N is the required number of modes, 2(N + 1) is the array sizt 
required. "NRAPH" calculates eigenvalue solutions by the Newton Raphson method. The forms ol 
equations solved are: J 

0= F-Ztan (Z) 

0 = F + Z cot (Z) 

where F is some complex constant and Z is the required solution. 


Subroutine "MRT" 


This routine accepts a filtered set of zero Mach number eigenvalue solutions and obtains solutions 
for the required Mach number. These solutions are then ordered according to least attenuated 
mode. For Mach number specified zero, the Mach £ 0 calculation section is skipped, and the only 
function of this routine is the ordering process. 

The Mach * 0 solutions are accomplished in incremental Mach number steps. The number of 
increments is based on a Mach number increment nearest .05 which subdivides the specified Mach 
number an integer number of times. The eigenvalue equations appropriate to the duct geometry and 
number of lined walls are solved for each incremental Mach number. The Newton Raphson method 
is employed with modifications necessary to restrain the solutions to the first quadrant of the 
complex plane. 


Subroutine "HWDBS" 


This routine calculates the eigenvalue solutions for lined wall-hard wall combination modes having 
hard wall mode numbers greater than zero. The soft wall solutions for zero hard wall mode number 
are used as starting points. The Newton Raphson method is employed for the solution. 
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Subroutine “RORDER" 


This routine accepts lined wall eigenvalue solutions in any order and orders them according to the 
largest real part of the downstream propagation constants. These constants are calculated in 


“RORDER". 


Naming Conventions for Transfer Variables 

Variables are transferred between routines through “common" statements and through the “call" 
statements. The naming conventions for these variables are given here for program reference. 

An unnamed common statement is used for the “mainline" routines and subroutines “INPUT", 
“ATNSPL", “IMPQD", and “IMPPD". In the “mainline" and "input" routines, all names are 
identical and are as follows: 


A AM, ADHS, 

POAC, DCON, 

DIAC, THKC, 

MTHKC Arrays of input variables of which each array member 

corresponds to a particular duct increment; the array names 
correspond to Mach number, lined wall duct height, perforated 
sheet fractional open area or polyimide number of plies, lining 
backing depth, hole diameter, material thickness, and boundary 
layer momentum thickness, respectively. The last three arrays 
pertain to perforated sheet linings only. 


OSPL Input sound pressure level spectrum. This array remains 

unchanged throughout the program. 

CFRQ Standard, 100Hz bandwidth center frequency array (200 Hz to 

10 KHz). 


SPL Attenuated sound pressure level spectrum. 

TSPL Three-dimensional array of modal sound pressure level spectra; 

the first index refers to the hard wall mode number, the second 
to the lined wall mode number, and the third to the frequency at 
which these modes are calculated. 


NSPL 

H 

RMS 

REORD 


Number of input sound pressure levels. 

Number of lined walls. 

Impedance option; RMS or discrete particle velocity treatment. 
Modal ordering option. 


35 



OVERMOD 

TL 

DL 

DHH 

G 

TT 

PT 

AN PLY 


Option to treat fixed number of modes irrespective of cut,, 
considerations. 

Total lined duct length (inches) 

Duct incremental lined length (inches) 

Hard wall duct height (inches) 

Thermodynamic ratio of specific heats. 

Total temperature, degrees Rankine. 

Total pressure, psia 

Number of polyimide plies or P e ^ or f ed e ^^f i0na ' ° P ' 
area if these parameters are constants, otherw.se zero. 

Number of lined increments. 

, 4tr> c ro” indicating the -lowest ai 

,ndex constants for Motion, will i 

highest frequencies for whicn eigenvd.u 

obtained. 

Static pressure, psia 

Velocity of sound, cm/sec 
2 

Density of air, gm/cm 

Characteristic impedance of duct environment (cgs rayls). 
Maximum number of lined wall acoustic modes to be treated. 
Maximum number of hard wall acoustic modes to be treated. 
Lined length of treatment completed at any point in the progro 
Soft wall duct height of first increment. 

Static temperature, degrees Rankine. 

.e only differences in naming the ^°T,?^Q S y a ^' 11 a e ^d "jDU^I^esprebvely^ ln^'ATNSPL" 
^r^of^uU and have no relationship to the va 

ed in the "MAINLINE" and "INPUT" routines. 

-iMpnn" and "IMPPD" differ from the above representation in 
ZSSSSSi ar e Sed, S all other variables are unused: 


LOW, IUP 

PS 

C 

RO 

PC 

NSWM 

NHWM 

X 

DHSI 
TSR 
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I 


"CFRG" (renamed "FQ") 

"SPL" 

"NSPL" 

Variable transfer through the "call" statements occur between the following routines: 


From "ATNSPL" 

to 

"GETV", "GETVPP", "IMPPD" "IMPQD", "UX", 



"RORDER", "HWDBS". 

From "UX" 

to 

"RTZ" and "MRT". 

From "RTZ" 

to 

"NRAPH" 


The variable names and meanings for these transfers are as follows (names in parentheses indicate 
name conventions in the called routine where differences exist). 

Calls from "ATNSPL" 

"GETV" and "IMPPD" 

POA (PLY) Number of polyimide plies 

V(VP) Particle velocity (cm/sec) 

AM Mach number 

RO Density (gm/cm/cm) 

TSR(TSDR) Static temperature (degrees R) 

D(BSP) Lined wall backing depth 

C Speed of sound 

PC(ROC) Environmental characteristic impedance (cgs Ray Is) 

FTPI Logical variable indicating first frequency pass through IMPPD 

FRQ(F) Frequency 

TSP(TSPL) Real variable indicating method of particle velocity treatment 
ZZ Calculated impedance 
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"GETVPP" and "IMPQD" 


The variable names and meanings are the same as for "GETV" and "IMPPD" wi: 
the following exceptions: 


POA 

TT 

PS 

THK, DIA 
MTHK 
"UX" 
FRQ(F) 

AM 

H 

C 

DBN 

MCUX 

Z(ZS) 

NN{NS) 

DHS(S) 

FSMAL 

ZZ(Z) 

"RORDER" 


Perforated sheet fractional open area 
Total temperature {degrees Rankine) 

Static pressure (psia) 

Perforated sheet thickness and hole diameter, respectively 
Boundary layer momentum thickness 

Frequency 
Mach number 
Number of lined walls 
Velocity of sound 
Modal attenuation rates 

Number of lined wall eigenvalue solutions required 
Modal lined wall eigenvalue solutions 

Variable indicating sequence of symmetric and antisymmetric mot:- 
Lined wall duct height 

Combined admittance, frequency and geometry parameter 
Lined wall impedance 


The variables are identical with those of ^LIX** immediately above except for 0 ■ 
following: 

AK cj/c 

G 


Z 


Duct height parameter 

Modal lined wall eigenvalue solutions 



"HWDBS" 


2(R) Lined wall eigenvalue solution 

DB Model attenuation rate 

N Variable indicating symmetric or antisymmetric mode 

AM Mach number 

FRQ Frequency 

C Velocity of sound 

K{M) Hard wall mode number 

DHS Lined wall duct height 

DHH Hard wall duct height 

FSMAL Combined admittance frequency and geometry parameter 

H Number of wails lined 

Calls from "UX" 

"RTZ" 

ALZRT Zero Mach number eigenvalue solutions 

MCUX(NALZ) Input is number of modes requested, output is number of eigenvalue 

solutions supplied 

{Variables "FSMAL" and "H" are identical with those of "HWDBS" above.) 

"MRT" 

ALZRT(ZRT) Zero Mach number eigenvalue solutions input as starting points for 
„ nonzero Mach numbed eigenvalue problem 

RK ulc 

All other parameters are as defined for inputs to "UX" and retain the same names 
except "MCUX" which is named "NROOT" in "MRT". 

' Calls from "RTZ" to "NR APH" 

FSMAL As defined for "HWDBS" 
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Z(ZS) Starting point for eigenvalue solution 

ALZRT(ZF) Eigenvalue solution 

The variable N' in routine "IMRAPH” indicates that the solution is required for 
antisymmetric mode. 


a symmetric or 
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APPENDIX III 

PROGRAM LISTING AND SAMPLE PROGRAM OUTPUT 
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r 3 
G 

>-) 

Q 2? 
o £ 

^3 


& 


w § 


c 

c 

c 

c 


***% NTISF $IJPPR FSSIHN - MUlTIMOnp ANALYSIS ~ 

SINGLE LINING POIYT Min HP PERF SHFpT 70MQ17 0 n k -7 
VAR I API F CHARACTER T STIC S ALONG DUCT LENGTH 

nt mfnjS l ON DP0t0O9) 


C 


c 


c 


rc al MTMK t MTHKC 

COMMON AAM(50|* ADMSISOI* POACHCI* OCQN(SO)* p 

i THK C I 50 » , MTHKC I SO) 

> * PSPM99), CERQI99), S p tt QQ )» TSPL(**?0*° 

* NS°L ♦ n, RMS* RPORDf OVP M 0D * Tl* 01* OMH, 

v TT, or, AN PI Y * J * LOW « TUP. PS. C. p 0. 

S NSWM* NHWM, X* OH SI* TSP 

HI MP MS I ON PLY(P^) 

FOUfVAt.FNCF (n°A N D * D 0A* P04CI1I* PLY(I)) 

1 t OHS, aohsi u ) * < D.nco 

* i.j ■» r 1 u r- «i|t # UT 


Pi 


( AM 




POUIVAl FNCF 
T 
? 


( 7 MACH, A A Ml 0 > ) t 
< vo ,nC 0 N I?) ) * 
C v M thk ,MTHKC( 2 ) ) 


IVOHS* A0HSC7J1* (VDPA NP* ° n 
CVOIA, otaci?) ) ♦ (VTMK . th 


CALI INPUT 
X = DL 
IOL = I 


Ohs 1 - dhs 
GO TO 40 
o s x * x ♦ ■ nt. 

IOL = TOL ♦ 1 
TP ( VM AC H .ME . p ) 

IP IVOHS .NP, 0.0) 

IF ( VO P A NP • N f . P. r > 

IF ( VO .NF . r . 0) 

ip rvniA .NF. r *. o) 

IP I VTHK .NF. r * A ) 

I C ( VMT HK .NF . 

jc I VM A T H .FQ. ^.^1 
PS = pt /( 1 .o + ( r * -1 


AM 

= 

A A M 

HDD 

OHS 

s 

ADHS 

unu 

POA 

3 

PQ AC 

( 1 nu 

T 

= 

OCHN 

lint) 

01 A 

= 

HI AC 

f l nu 

THK 

- 

THKC 

unu 

M TH< 


MTHKC (TOL) 


GO TO 40 

.<>) *AM*AM/? .0) **IG /(G -UD) 


TS P = TT/U.'VIG - 1 .ni*AM*AM/?.D) 


TS - PT?.l S+S.9/o. n *I TSR-A91 .69) 

C * ' J 314S.0*S0RT< TS/?^3.16) 

RO = 27^.1 6/T$*PS.M .01 3?5E6*T .293E-3* £R Q 47. 


pc = r n*c 

40 CALt ATNSPL 

00 SO 1 = L 0 W * I IJP 

so noo( n = ospl(I) - split) 

WRfTF I**F) X 

WRITE U, 7 ) (C FPOI I T.OOOI n * I =LOW, I UP) 
WRITE M , F> X 

WR HE I I.T) KFRQl n *0301 n .! =LHW* ? UP) 
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|F (ABS(X-TL) .GT. SI) GO TO 38 
WRfTF u f m 
write n,i 7 > 

FT&ST F = !O r *Lnw * 100 
fFIOST = LOW 

80 LAST F = MIN0(FT9ST F * *00, 100MUP ♦ ICC) 

HAST = MINK TF T° ST * 3, I U° ) 

W3 T T F (fc,R11 F 1 ° S T F, LAST F, (SPLIT)*! = I F T R.ST , ! L AS T ) 

WR1TF (l t P|) FIRST F, LAST F v I SPLIT >,I = IF! PST, HAST) 
fF ( IL AST . EQ. I UP 1 GO TO 90 
FfRCT F = LAST F ♦ ICO 
TFJRST = ILAFT ♦ 1 
GO TO 8C 
°0 IM ) = I UP -1 

01 70 T = LOW, f Ml 

oro( n = coRom + oantm))/?* 
to cffo(I) = icfpqiti ♦ cfrquh u /2. 

WRfTF ( If IS) X 

WMTF I 1, l 7 ) (080(1 1 , I =L OWf I Ml ) 

STOO 

5 PI rm AT (*? 08 DIFFFPFNZF AF TFP • , F6, 2 ,* TNCHFS OF I INING 1 ,//) 

7 FORMAT ( • SFR.O.fio^j 
11 FORMAT (M FRPOUFNCTFS FINAL SPECTRUM*) 

IS FORMAT (*0 C R 0TFFPQFN3ES (SMOOTHED) AFT FP * » Fft . 2 , * INCHES (IF* ♦ 
l • LININGS//) 

1 7 FORMAT ( * * , SF1C.2 ) 

81 FORMAT ITS T A , * tq* t 16, AF11.2) 

FND 
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c 


SUPPO'jTTNF TNPUT 

RE At MTHK, MTHKC 

INTEGER FIRST F 


), OCDN(ftO), OIACt 


CGMM ° N AAM(ftO), AOHS(5:>), POACf 

TM<: (SP| , MTHKC(SO) 

» osplcoo), cfrocooi, sni^i, tspi (6, zo, oo> f 

NSPL, H, RMS, RFGRD, CVPMOD, TL, 01, OHM, G, 
TT, PT f ANDLY, J , LOW, IMP, PS, t , RO, PC, 
NSW^, NHWM, x t OH SI , TSP 
OIMFNSION PLY(ftP), SPL NAM(JP) 

F3U f V AL FNC F <0™A MO , POACH), PLY(l)) f (AM, A AM ( ] ) ) , 

1 ( OHS, AOHS( 1 ) ) , CD, OCONH ) ) , (HI A, DIAffl)), 

f THK, R , THK: ( 1) | , I MTHK, MTHKC f J ) ) 

COMMON /NJN/ NJNI^fOg) 


C 


C 


FQlJIVAt.FNCF (VMACH, AAM(2)), (VDHS, 

1 <V0 ,0C1N(?U, ( VD I A , 

? (V M THK , m THKC(2) I 

SFT CONSTANT VALNF F L AT, S ( V A R I A T I ON 

VW ACH = "*o 


AOHS (? ) ) , ( VOPA NP 
OTACC?n, ( V THK 

AtONC OUCT LENGTH 


VOHS ^ o*o 

VO p A NP = 0*0 

VO =^*0 

VO! A = C • ° 

VTHK - * m n 

VMTHK = 0.0 


poAcr 

THKC ( ' 

rF NQf, 



S3 


fes f—P 

Pis 
v © 

■J fM 

'3 as 


P: 

f^S 


o 


OIMFNS IPN T I T(_ F | 1 R | 

RFAO (5, A ) T IT?. F 

(S * ?3,4 * SWM * HWM * H, PSPL, RMS, "FORD, nvRMOO 
NSPL = A 

TF (PSPI . FQ . O.0J r,n TO 44 
no 43 f 3 I, NSP| 

43 SP|( I> = PSPL 
GO To 4 5 

44 PFAO (ft, 4) S°L N&M 

READ (5, 23) ( S p L ( I ) , I = 1, NSPL) 

45 IF (HWM ♦ F Q • 0*0) Uf^M as ],0 
HWM S AMfNl ( ft . o , HWMJ 
NHWM = HWM 

IP (SWM .FQ. 0,0) SWM = 1.0 
SWM = A M INl(?0.o, SWM) 

NS WM 3 SWM 
K = J 

00 4 ft J = t , ?o 
DO 46 T = 1 , N $ PI 
4ft TS Pf (K , J, I ) = r *n 
CPRQ(l) 3 ?oo. 
on ft ft I x 2 ♦ NSPL 
ftft CFRO(I) 3 c FR Q ( I- 1 ) ♦ 100 * 

RFAO t ft , 2 7 ) Tt , 01, OHM, G, TT, PT, A, p 
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o u o 


j * (tl ♦ c.n /pl 

READ ANPLY, AM, DH S « OPA NP t Pt OfA, P, MTHK 

WR ! T F t ft , 1) T [ T{. F 

WRITE (1, ft) TITLP 

IP ( AN P l Y . FO . C. n ) WRITE (ft, 2) 

TP (ANPIY .ME, P.P) WRITF (ft, 3) 

WRITE (ft,*') 

mw = A/icr. -.9 

IF (LOW .FQ. P) LOW = I 

RIP = B/IOO. -.P 

IF ( (UP • FQ . o ) IIP = NSPL 

WRITE (ft, 9) N$PL, NSWM, NHWH, H, PSPL, RMS, PEOPP, DVRmOO 

WRITE ( 1 , P) N$Pl, NSWM, NHWM, H, P$P1 , PMS, RFOPP, OVRMOP 

WPfTE (6, *) TL, PL, THH, G, TT, PT, CFRG(IOW), CFPQ(HJP) 

WRITE U, P) TL, PL, PWH, G, TT, PT, CFPQILOW), CFPOCItJP) 

WR(TF (A, 1C) ANPLY, AM, OH S 
WRITE (1,10) ANPLY, AM, OHS 

IE {ANPIY .N c . P.o) GO TO ftfl 

WR ITE (ft, n ) OPA NP 

WRITE ( 1, II I OPA NP 

WRITE (ft, 14) D, PfA, R, MTHK 

WRITE (1, 14) P, PI A, R , MTHK 

GO TO 69 


63 WR ! TF (ft, 12) OPA N 0 » D 
WR ITE (I, 1?) OPA NP, P 


69 

FJFCT * AM IN 1 

( AB S ( AM) , 


OHS 

, 

OPA 

NP, 

0, DIA, 

Pt 

MTHK ) 


IF (AM .NF. 0. 

0 ) 







GO 

TO 

70 


REAP (5, 27) 

( AAM 

( I ) 

, 

I 

* 

It 

J) 





WRITE (ft, 30) 

( AAM 

n ) 

t 

I 

= 

1 , 

J) 





WRITE (1, 30) 

( AAM 

U 1 

, 

I 

s 

It 

J) 




70 

TF (OHS.NF. 0. 

) 







GO 

TO 

71 


PFAP (F, 27) 

( APHS 

( i ) 

, 

I 

= 

It 

J) 





WRITE (ft, 3U 

( AOHS 

( ! ) 

, 

I 


It 

J) 





WRITE (1, 31 ) 

( AOHS 

( T ) 

♦ 

I 

= 

1 t 

J) 




71 

IF (ORA NP .NF 

. O.o 

) 






GO 

TO 

72 


REAP (5, 2 3) 

( POAT 

( T ) 

♦ 

I 

= 

It 

J) 





IF (ANPIY .NF. 

P.C) 







GO 

TO 

77 


PERFORATED 

S HF F T 

’PER 

CFNT 

’• OPEN 

AREAS 




INTERNALLY, f p AC T I ONA L 
FXTERNAILY, RTR CFNT 
WRITE (ft, (ROAC (II, 1 s 1, J) 

WRITE (1, 3?) (POAG (T), l = 1, J) 

GO TO 72 

C POLYIMTD - NUMRER OF PLYS 

77 WRITE (ft, 37) (ply (I), 1 = 1, J) 

WR I TF ( 1 , 37) (PLY ( I) , ! = 1 , J) 

77 IF (D .NF . 0.0) GO TO 7i 

READ ( c , 23) ( PC ON (I), I * 1 , J ) 

WRITE (ft, 23) (nr ON (II, I = 1, J) 

WRITE (1, ? 3 ) (PCON (I), ! = I, J) 

73 IF ( 01 A .NF. P.o .PR. ANPLY , NE . C.C) CQ TO 74 




R F AD ( 5* 77) ( D T A C (I), I * It J» 

WMTF (6, 34) (DTAG (Tit ! - l, J) 

WHTF (1* Ml I (MAC IT)* T = 1* J) 

74 1 F { R .NF. 0.0 .DR. AMPLY . NF - O.f) GO TO 75 

P F AD (5* 2?) I TM^C (Tit I = 1 * J) 

WRTTF (6, 3 s ) ITHKG f T I * I = l* J) 

wp t tp n * 75 » i thyt i r ) » f - !• J) 

75 IF IMTHK .NF, n • * .OR. ANPI Y • NF • 0.0) GO TD 76 

PFAD (5, ?3I (MT HKnil, I = 1* J) 

WRITE (6, ’5) (MTHKC.U)* T = 1 t J) 

WPJTF ( 1 T 36) (MTMKCtn* I = l* J) 

76 l r (PSPt. .NF. 0.0) GO T n 130 
IF f EJECT .NF. r.'M GD TO 7P 
W» I TP (6* 6 ) 

WRTTF 11*6) 

7q W« I TF (6*17) SPI NAM 
WRITE (1* 17) SPI. NAM 
FIRST F = 700 
! F IP ST = 1 

PO l A ST F * mjn'M^PST r + 700* 100*NSPL ♦ TOO) 

(LAST = MfNDUFTRST *■ 3* N$P1 > 

W3(TF (6, PI) FIRST f 9 LAST F, (SPL(I),! = (FIRST* (LAST) 

IF (LOW .LF. (FIRST .AND. (FIRST .LF. !U° .OP. 

1 LOW .LF. (LAST .AND. J L A S T .LF. (IIP) 

7 WPITF (l, pi) FIRST F, LAST F, (SPLUI* I = I F IRS T f (LAST 
IF ( (LAST .FO. NSPL ) GO TO l?'* 

FIRST F = LAST F ♦ ICO 
IF IP ST = (LAST ♦ 1 
GO TO PC 

1 70 WR TTE (6, 6 ) 

WRfTF (1,6) 

130 CONTINUF 

PS= PT /(l.C+(G -1 ,0)*AM*AM/7 .0) **(G /( G -1.0)1 

T$R = TT/C1.0 ♦ IG - 1 .0)+AM*AM/?.O) 

TS * 277.1 6 ♦ 5 .0 /9.o*( TSP - 401.69) 

C = 73(45. DTSOP T( IS/273, 16) 

PO = 273.1 6/TS*PS/l .Cl 325E6*1 .293F-077 />RP47. 

PT = R 0*C 

WRITF (6*19) 

no FO I s LOW* TUP 

SW =(?.?*CFRO( TI^DHS* 2.54/C 4 1.0) / SOPTd.O - AM * AM ) 

MW = ( 7 . *CFP Q ( H *DHM*? .5 4 /C ♦ 1.0) / SOR T ( 1 . 0 - AM* AM) 

IF ( SW .GT. SW« .^R. OVRMOD .GT. 0.0) SW * SWM 

IF (MW .GT, HWM .OR. OVRMOn .GT. 0.0) HW = HWM 

TH = HW 
IS = SW 

IF ( IH .(T. II TH = 1 

IF (IS .IT . II IS = 1 

A a IS 

A = S^LII) - 1C.*4L0G1D( 4*1 H) 


...xOJJ 1 ® 



NJNC 1* 1 J - tS 
NJNI2.I) = 
nn 2 ? L « If IS 
03 27 K - t*I M 


?? TSPLtK til n = 4 

write 1 6*?i. ) cfrqi n » 
so nspt m s 5f>im 

RFJURN 

1 FIRM M ( • 1 *$*S *F * * 1 BA 4 I 

2 FORMAT <«C R FC TANGUL A* 
1 FORMAT ( *r Rrf TANGULAR 




AN&1. YSI S 

analyst s 


NJNl? *1 * * T$PL< 1 tl * H 


» ppofORATFO SHEET MOPFL') 
- POLYfMTOF MODEL') 


4 FORMAT 

5 format 

6 FORMAT 
q FORMAT 

SF14.6* 
FF 1 4 *6* 
RF1.4 .6, 
9F14.4* 
7F1 4, 6* 
C.F 14 .6, 
OF 1 4 • 4 * 
q FORMAT 
G 114* 
M 1 1 4 * 


• s* // 

MOOES 
MOPE S 


( 1PA41 

f«0 TN p U T PARAMETERS') 

IM f .lRA4J 

(.C<, FU.».< 6X,*TDTAL DUCT LFNGTH* / 

6X, • TNCP FMFNT LENGTH'// 

6X,'HARP WALL DUCT HEIGHT'// 

6X, 'GAS CONSTANT » / 
fX, 'TOTAL TFMPcq ATUPE • / 
feX* ' TOTAL PRESSURE'// 
fX* M OWE ST FRFDUENCY rONSIPFPEP'/ 

f X f ' H TGHFST FREQUENCY CONS! PE RFC * / ) 

♦ ft X » ' NUMB FR OF S p L 

fex * • m a x I m um number of soft wali 

f^X t 'MAXIMUM NUMOFR C3F HARO WALL 
, F'4.6, fcX,'MU^E« OF WAILS LINED*// 

JFLA.S. 6X,*UNIF0»w SPL •// 

6X,«BMS PAOTISLE VELOCITY*/ 

FXf 'REORDER OPTION'/ 

6X*'0VFPRT0F MODF OPTION'/) 

F 1 3 ♦ 6 » 6X, 'MODEL SE l E F TOR * / / 

6X, 'MACH NUMBER ' // 

6X,'S0FT WALL DUCT HEIGHT'/) 

(RPF14.6* 6X» 'PFRCENT OPEN AREA'/) 
f F 1 4 « 6 * 6X f 'NUMRFR OF POLYIMTO p l y S'// 
6X, 'RACKING SRACE'/'2M 
(IX, 1 R A 4 / FREQUENCIES 

( F! 4. S< 4X, 'RACKING SPATE'// 
ftXt'HOLE OTA ME TER' / 

6X, 'MATERIAL THICKNESS'/ 

6X* 'MOMENTUM TH l CKNF S S' / * 1 ' I 

*»p foual PR fssure MODAL SP1 
»p frfo number o c number 

SO r T WALL HARP WALL 
{IX* FS.o, |6* III* FIR. II 
{ RF R • 1 1 

( «o MACH NUMBERS 

SOFT WALL DUCT HEIGHTS 
P FR C FNT PPr N AREA 
RACKING SP A C F S 
0 I AMFTFR s 


CONS I PERT D* / 

C ON SIDFRFP* // 


KF1 4. 6* 
LF14.6, 
MF\4. bt 

10 format 
l FI 4 * 6 » 
FF14.4, 
u FORMAT 

12 FORMAT 
1 F14.fc* 

13 FORMAT 
1 4 FORMAT 

Y FI 4 • 4* 
AEl 4 , 6* 
ZF| 4,6* 
19 FORMAT 
\ 

2 ' 

21 FORMAT 
23 FORMAT 

30 FORMAT 

31 FOR m AT 
>2 FOR M AT 

33 format 

X4 FORMAT 


SOUND PPESSHRF LEVFLS*) 


' / 

PF 


SPL FOR' / 


EACH MODE ' /4 X,R ffeK ** MODES* )//> 


( 'P 
( 'P 
{ *P 
( '0 


« / |9F 1 2.4) ) 

• / (9F1 2- 4) ) 

* / (2PF 12 .4 I I 
•/(RF12.4) ) 

' / {9F12.4) I 
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1 5 
P 6 
77 
P 1 


PQRM AT 
FOP* AT 

FORM AT 

F0R m AT 
fm n 


«/|QFl?« 6 l) 

THTCKNFSSFF kirrccc */IOFl?*6>l 

; ^'"^or b vn, n pl y S . / « «■=« ■ • 

\,r,, 14, • TO', 16. 
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c 


c 

c 

c 

c 


SUBPOUTINf ATMSPL 

MTHK, MTHKC, PI/3. 14159/ 

INTEGER FT PST p 

COMMON A A M ( 50 ) * A0HST5P), PnAC<50), DCON(50), OTAC<50), 

l THKC (50 ) , MTHKf ( 50 ) 

7 * nSPJ.(9P), rF90(99), St>L(Qn) f T 5 PL ( 6 , 2? , 99 ) , 

1 NS°l, M, RMS, PEORO, OVRMOO, T L * 01, DHH , AKOUM, 

4 TT * PT t AMPLY* JDUM, LOW, (UP, PS, C, PO, Pf, 

5 N S WM , • NH WM , x, OH SI , T$f? 

OTMFNS TON PLY( 5 r > 

EQUIVALENCE f HP A NO, p 0 A , POACH), PLY(l)) f (AM, AAM(U), 

1 (OHS, A OHS ( 1 ) ) , (0, nrONCl )) , CPI A , 01 AC { 1 ) I , 

7 - t THK, p , THKCdl ) ♦ (MTHK, MTHKCdH 

EQUIVALENCE (VMAfH, AAM(?|| V (V0H$, AOHS(?)), (VOPA NP , POAC(?))» 
1 (VD ,DCnN(2U, ( VD I A , Of AC (?) I , (VTHK , TWKC(M), 

? ( vmthk, MTHKC(2I ) 

COMMON /NJN / NJN(2,oq| 


COMP! EX ?A A?) ,F$mal v FS, 7Z, SAVE 17 
D [ MENS T CN D AN ( 4 ? > , N N ( 4 2 > 


V = 40 * 


FTPf = .TPU C . 

IF ( RMS .L E . o.) CO TO 50 
WRfTE (6, 19) 

19 FORMAT (MM 

fP (AMPLY, NF • r ‘,°)fAl. L 3FTV(P0A,V 
CALL GFTVPP (POA, v, AM, TT, PS, 
50 00 70 I = LOW, IMP 
TSP = SPL ( I ) 

IF (RMS .IF. O. ) TSP = -TSP 
FRO * CFROf I) 

TF (ANPIY * NF . O. r ) CALL [MPPO 

1 

CALL TMPQD 

1 TT, PS, THK, 

55 WR !TE( 6,9) FPQ, 77 , V 


AM, RO, TSP ,D ,C ,PC , 4-50) 
THK, P I A , MTHK, 0, f, PC) 


( F TP I , FRO , POA, V, AM, RO, 
TSP, 0, C, PC , TSP, 77 , 4-55 ) 
(FTP), FPC , POA, V, AM, 

OTA, MTHK, 0 , c, PC, tsp, n I 


A = 0.0 

NS =( 2 . ?*FPO*HHST* R ,54/C ♦ 1.0) / SOP T ( 1 * o - am*am) 

NH "( ? . *ERQ*PHH*? .54/: ♦ 1.0)/ SQRTtKO - AM*AM) 

TF (NS .G T . NSW* # OR . OVRMOO . GT . C.O) NS = NSWM 

TF (NH • C»T * NHW« .OR. OVRMOO . G T . 0.0) NH = NHWM 

IF (NS .LT. 1) NS = 1 

IF (NH .LT. 1) NH = 1 

MCtlX = NS 

CAM UX (FPQ, AM, m, 0 BN , MfUX, Z, NN, OHS, F S M A L , 17 ) 
WRITE (6,5) ( 7 { J),09N( J I ,NN( J) , JM ,MCUX) 

A< = 2,*PT*CFP0( I )/C 


REPRODUCIBILITY OF THE 
ORIGINAL PAGE IS POOR 


4-1 



!?.£ -•» CAU RWOE* IA*,AK,G.»rUX,7.NN.OBN> 

H = MINO( MCUX, N S, MJNll.m 
NH = MthTlNH, NJNf2tlD 

AlMAGtm ” - NE ‘ A,MaG(n? ' 

\ GO TO 21 

W3fTE ( *t ?2) CFRO(T) 

OftNtll = c . v -t pin 1. *♦ EPPOR. DUPLICATE SOLUTIP 1 

?? 1 ” PM&T nh" ATTFNUATJOn'aPPlIfD 1 FOR fW OR TWO “OOFS 

?1 WRITE { 

mfu;n'= 

IP (NH .EO. 1 1 r ’° Tr * 

^fE < 0 RNfJ! .^ n - G0 T0 ? ° 

■ ^ 0 HWn^ ? ;n J ,.0,,N.A M .FRO.C.K,OHS,OHM.FSMA L .H.FS, 

WMTF Ifc.T) J.K.ZIJI.03.N 
10 ts°l ( k » j * t i - Tsn(«.J.n - 0B * 0L 
?0 CONTINUE 

NJ5 r M IMP(NS # NJNUi!)) 

OT “Jf! J = liNS 

DD BO K - t *NH . , t 

•>0 A * A ♦ lO.**»TSOt (K,J,n/lP.» 

' cpi m = ir.MLnr.iotA) 

IF (MOOU-LOW, 71 .F3. 1’ W * ,TF ,6 ’ 

td continue 

WP ! T E ( E»n * 

FIRST E = PCD 

*0 LAST^ =’mNAIFI^T F ♦ 5 ^, 100 *NS»l. ♦ 1^1 

^’iSPLUI.T - IFIRST. UAST, 

IF ( IL A ST .FO. NS°L 1 RFTIION 
FTO ST F * LAST F «• l 0 ' 1 
lFt°ST = T L A S T ♦ 1 

r,n th rd . . 1 ~ . 

p ! ^ Jn'^FOUFNC^s' * nfw spectrum aftff., fa.o, 

\ • INF HF S OF LINING*! 

? prjpM AT ( 1 p M 
^ FOP MAT | «0« ,1PRFl*.*> 

■> format (*F onnTS an" ' DRN nn* i ' 

^ » t 70 P AL * * MAb 

7 2I1PF1A.F, MO.S, OPFR.F, I3M 0 RN NN* 1 

,, rnOUHT EW M HWM 7REAL 

\ f F ORMA T T IlF. IA. 1 - 14 . 5 . Fn. 5 , PPFF.A, HI 
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9 FIRM AT < f C 
1 * 
FNO 


t /% FpfQHFNCY Ffc.Ot * H7’/ ZR*ZI 
v* « F14*51 


f lP?E1^.5f 
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Simp OUT INF IMP orjprpf t F , OLY f VP f am, ro.tsdr, rsp,c* roc f tspi , zz , », 
C IHPFnANCF - PHIYIMITF ZDMQp 

COMPLEX fMOLX f /Z 
LOGICAL *4 FTP T 
DATA P T /? * l 41 *<5/ 

COMMON RUN CHI 4 4 9 ) « F0(99) # SPl (99) , Mf)RF III R R o j t NSPL, T HE , R F<; r r 
niMFNSICN PSPL(PO) 

COT ( X)=fC!S(Xl/SINf X) 

D* R S P 

f ABSf T<;c>L 1/20.) 

TKD= ?.* P l *F /C *q*? . S4 
I IN = ^ 

IN = 0 

IF( .NOT .FTP T) GO TO Rn 

u=3. l7E-P/?.ooq* T SOR**l ,5*73 4.7/(TSOR4-2U. ) 

RVnc = 464C.*U** > L Y**l .333 

V® GF*. C5*r* AM* AM /(R VOC /ROC*- SQRTf ( R VO C/FOC ) **? ♦ . 7 1 *PLY* P LY *AM*&v 
IF ( T S p L .LT. o # *| WRfTF (f,l) VPGF 
\ format < »o GRAZING Flow PARTICLE VELOCITY = *, FU.4 f * FROM SUP°- 

I •niJTTNF I m d of) * ) 

FTP! = .FALSF. 

50 CONTJNUF 

p ! AC = 2. Ff.*RO *PL Y*PL Y*VP 
PVAC*U*.?0*5*pi v*d|_y*f 
RZ *P f AC*P VAC +° vor 

RX*R0*F*t. e A9*PLY** 1. 5-5. 12F-R*PiY**3*vP) 

I F ( T SPL • CF . 0 • I GO TO *0 
R 7 P V = R I AC/VP 

PXPV=- 5. 1 2F- 5*R p*°l Y**R*F 
IT =TMPL X(PZ ,P X-POC*COT( TKO) ) 

OAR = CARS! 77) 

SOP TV = SORTfVPGF**? «. (P/OAR)**?) 

FNR * V° - SOP T V 

,PPQ7 = -2. /QAR**4A( R7*p jpv 4- R X*R XPV ) *P* P 

FNRP = 1.0 - PP07* . 5/SQR TV 

VPN* VP- FNR /FNPP 

IF( VPN.LT.O. } GO TO 4 I 

IF ( ( 4RS ( VPN-V° ) P0C1* VP| . I T. o # ) GO TO 60 
I TN = ITN + 1 

IF ( rTN.GT.50) GO TO 7r> 

VP * VPN 
GO TO 50 
41 VP=1.0 
' IN = I N+ 1 

IF( IN.GT.3) GO TO 70 
GO TO 50 
70 CONTINUE 
77 = 7 7 /ROC 
RETURN 1 
60 CONTINUF 

II * 7Z /ROC 
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RO 


5 


17 


VP*VPN 

R E TURN 1 

T<D=? # *Pl*F/f*RFf>*?.«;4 

zz=cmplx(p; /» oc t^x/pnr-cnrf tkd) i 
RETURN 1 

FNTBy RF T y (PI V, VP, 4«, »n, TSOR, RSP, C, ROC, *) 

ST * T = 1 , NSPt * 

p< ^ p lU> = 2. P£-4*i*.** ( ( I ) /20. • 

nru wl7E ' e/;> *'‘ Pq,5tTSr)R * <tl *5*734.7/ (TSPR+21 6.) 

DCM s BSP * ?,5A 

RVDf=464^.fc{j*PLY*if'l .^^3 

wRt’Trnr!rvPGF /,pvDr/,oc+sQpTuRvDc ' pnr, *‘ ? ^ n * oLY * piY * flM * AM '> 

TK = 2*C*Pf/C 
VP = 

p 14C = 2. P6*Rf|*PL Y*Pf y* VP 
PP0?=0. 


15 


p S7= 0. 

m 15 I - 1 ? NSPL 
TK 0= TK $ FQ ( r J*DCM 
«VAC=U*.?^5*PL Y*oi.Y*FQ( I ) 
ft Z =R ! AC «-P V4C+R VDt 
PX-RO*FO( I }*( »54P*D| Y**x,5-*^ 
71 =fHPL *{R 7 , 0X-Rnr*CnT{ TKD) ) 
0AB=CARS( 77 ) 


12F-5+PL Y**3* VP ) 


P17F={ R$Pt { ! )/ OAR)**? 

PQZ = PQ7 •fPQ7F 

PZ f *-5 • ! 2F- 5*Rn*t>LY****FQn J 
PZR = 2. 8f«pn*PL Y*Pl Y 


PP02 =PPOZ-2./OAe**4*<*7*P7R+RX*PZH *PSP1 f T 
S^RTV=SORTC V p nF*VPGF+Pn Z) 

FA 1 — VP -SORTV 


*PS p L i I) 


F A 1 P* l * - PPOZ*. 5/SQR TV 

VP l =V p -FA1/FA1P 

IF( APSf V p l -VP ) . L T* .0001* VP 

VP = VP! 


PFTURN I 



* tt oc THK ■ 0 t A • HK f RS P * ^ * 

SUBPOUT WE IMPOOIFTPI.F.oOA.VP.AN.TT, S. 

{m^pIncf’-^pf'FOpateo s^eft , , mc « f «ii and • nspl, thf.bfst (2?. 

COMMON RUN'CH(A40), SVLI 

IMMFNSinN °s°l i Cf> > 

CO**PLFX CMPLX*77 
LOG IC Al ** FTOI 
PFAL MTHK 
DATA <Mn.lM59/ 
nrl X)=COS( X) /StN< x) 

p s? # 2£.. 4*1 r>.** ( ANSI TS° L ) / 2G * ) 

T < D l 7# *PI*F/r*RSO*?.F4 

!TN=n 

|N| =f> 

JFt .NTT .FTP?) GP TO 5P 

th=tt / 5 1 ^ • 

pp? = PI/?./SQB T t?*) 

2 F V“o?s1mi /Pn*..2-n/-o.**.l/E**«-r T2 *™ K,OU ’ 

S;o?.!’T7tTHi/OFL/Pn»*TH.TH/«TM*.M6) 

:^^^^.p^.g I& /-hk,*c.^-hmpv ? .sopt«pv ? - ? m.: 2 -.o,-ou 

^TMK )*F**VM* AM ) ) 

-• "" 

I « 0 U T T N F T m PQ3M 

pyp| = «FAL SF * 

so mntinuf^ t nlA * p * PH&/VP 

::^!p 3F--TH«"^S ^^MMTHK/niPM.-POP./IPOPPSOPT.OEf.TH 

1 4- 4 1 6 t M 

si s:i^wi5S ,».,-i.»-'w<> 

IPIT^PI cf.c.) go to ro 

, ’ 7 pv "■ if«/VP‘tl.*3.B*SN*SHl*«nC 

p*pV=-F X/C THK+Ol T )* 1 . PH*-*’ f ■ 

T e (olt .fo . ^* r > B yov * JJ; 

7 7 =cmplx( p ? f PX-pnr«r.nT( tko> i 
SORFV = Ae SORT(VPr.F«*2 ♦ ( P/OAB 1**21 

ppn7 VP /-°2 P ./o^*^*<P7*^P^PX* R xPV.^*P 

p M pp = i.p - ppnz*.F/s^RTV 

VPN=VP-FNP /FNRP 

i:!ri5 1 c^:i-“."0«tIvp|.iT.0.i on rn so 
itn=itn*i 


beprodugbility of raB 
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V D SF 4- 1.0 


I p ( ! TN • GT . F ^ ) GO TO 70 
VP *VPN 
GO TO ^ 

*1 VP=1.0 
f M = T N > 1 

IF (AH .NO. 0.0) VP = 

IF< IN.GT.O ) GO TO ?n 
GO TO *5 0 

70 CONTINUE 
W9ITF 16, 71) IN, I TNI 

71 FORMAT (» PET IRN C Q VIA 70 F«QH T MPQD IN = , ,T7, • I TN =•. 1 1 ) 

17 = ZZ/POC 

RETURN 
60 CONTINUE 
77 = 17 / ROC 
yp =vpn 
RF TURN 

80 TK0=?.*P!*p/C*9Sp*2.F4 

Z7=rHPLX(PZ/P0C ,P X/ROf.-COTC TKD) I 
RETURN 

ENTRY GFTVPP ( POA , VP, AH, TT, PS, THY, D1 A , HTHK, BSP* C, ROC) 
DO *5 I = 1 , N?PL 

S PSPU ! ) = 2.?F-4*!0.*M$PL( I ) /20 • ) 

REAL K = .OF ♦ . 1 1 * 01 A /MTHK 

TH=TT/F!S. 

P p 2 = P I / SQP T( P . 0 ) 

DEL = PS / 14 • 7 

E-l .0251*1 1 . /P0A**?-1 ) /PDA**.l /FXP( . FC72*THK/0I A ) 

rvoc = Rvnr/ippor) 

R V DC = • 0 77# THK/DEL*TH* TH /( TH*. 41 6 ) /POA /? . 0 

VPGE=0. F*P FAL K *C * 4 m* * M / ( P VDC ♦ SO RT ( R VDC **? + . 5 *RE At k*F*AM*AH) ) 

2 * O.F 

WRITF (6, II VPGF 
VP = 10^.0 
IN = 0 
1 7 PP OZ =0 . 

PDZ=0. 

IN = IN ♦ I 
DO IF ! = 1, MS PL 
TKD»?.*P!*FOU ) /C**SP*2 .*4 

RVAC=.rCCC7B3*TH**. 7 s*sqr T(EOn ) ) / POA/SOPT (DEL*< TH+.416 ) I 
1 * ( THX /HI A*1 .rt-OQA ) 

SN = 2 • C 4*TT A*F0l! ) *POA /VP 
R I AC*PP?*VP /C*F*E X° C- 1 . R* SN*SN) 
r; =r I AC ♦p VAC +R V nr . +p VDC 

DLT = AMAXM.R c ;*OTA*{l* n -.7*SQRT(P0A)| — .OPOO!2*VP/POA, r '.0) 

RX = •CCCAfcO*FQ( ) )* (THK4-DLT) 

RX = .CC0460*FOn )* <THK*Ot T) /POA /SQR T(TH) 

17 «CHPLX(R Z ,RX-rnT{ TKDM*ROC 
OA R -C A F SI Z Z ) 

PO?F^(PSPl I I ) ZQAR)**2 


.*p5 



Pnz=Piz*pnzF 

RXPV = -PX/{ THK+DLT)* .00001 2/POA 
P7PV=RTAC/VP*{ 1 . A + ?.6*SN*$N) 
py = PX - rnrc tko \ 

15 PPD7=PPP7-;?./OAP*fr4*(R7*R7PV«“PX*RXPV)*PSPL{T)*P < >PL(I)#Rnf 1 ! [ *? 

$3RTV=SQRT( VP0P*VP0F-*-POZI 

FA 1 = VP - SOP TV 

FA l P= 1 oo 0 7* . 5 /SOP TV 

VP 1 =VP -PA1/PAIP 

IP ( VP 1 .IT. 0.01 0,0 TO ?C 

IP(ABS(VPl -VP ) .1 T. . 000 1 * VP | RETURN 
VP = VP] 

GO TO ] 7 

?0 IF (IN .EO. ?} GT TO ?] 

VP = 1.0 

IP (AM .Nf. 0.0) VP = V D GF ♦ 1.0 
GO TO 17 
?l WRITE (6,??) 

??. FORMAT (* R FTI |P MFD VIA 21 FROM GETVPP IN = «. T71 
return * 

FNO 
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SURPflUT TfS/F IjXfF, am u 

COMP I FX*p An t > > * ** * 0BN f MCUX.. 7 c 

PPAL *4 PRM ( 4 ?) ’ f ^ SMai ,ALZRT ^?) * * 

ntMFM^ TO^ N ${ 4 ? I 

**** pk r S K= H /r 

P I*0.l41*c 

41 = CMP{_ * ( o r i r * 

0 HS x <; 

4 D=l ./7 

P ^ S«^L F r , ^t*4;i p T •' :l • 0, rif1 TO 50 
WP ITE (6,1 | fUa ?R ( T ’ MCU)f »H) 

" i.EAi,%Wr % v;; 1 " 

N = MClJX - j * GO 

" n * J«] t N 
l=Jfl 

00 4 k = i t MCUX 
4 R $ R F 
4BS f M - , 

JF ( ftflSP F 

* or . j . f- b \ qh to a 

= CMP <-*M.P?r f n.n, 


NS 


* St FSMAL, l) 


> J = 1 t 
TO 64 


*cuyj 


t-l fMCUX 

PF .GT. ,.f!,, 7 ^ , ^'; 4 ' M *G(A.7 0 T,K)), 
[ M *GT. t.F-5) no _ 


‘FS>«4l,U2PT,R K ,H.MCUX,OBN,AN,OHS,ZS,NSI 

1 A ^ t A ^ 


f p f ABS f m 
al;pti ji 
GO TO ? 

4 conti nuf 
3 CONTfNUF 
64 CALI MPT 
&F T(JRnj 

.m 

* •lipfi-NC, JsViri T ?f 9FAL ° ART THE ABOVE . 

^ * IJL AT f ONS A P F B YD A S ^ Fn cL^'^ NFGATfVF, $OMF r k\ r 

WPITf Cl, n?) SU in S A,F r ' IVAL,0 * , //25('I S * FF f°' ,F>,rY 4N0 p °IMTEI 

1 MMOFOAMrr', , . / ''' 2 F l r J' r ^ R F»f. PART OF THE • 

: ,,t4 T'ONS ARE BvPA^Fn FOR r ™^ .o^ AJIV ^ «M F 0410- 

c ULTS ARP f NVAL TO. *//?«%(, * , -OUENCY AMD PPfNTFD 


MCUX- = 
ffturm 
1 FORMAT 
1 

FND 


(,P 40M f TTA^CF °ARAMFtcd 

•0 PPELCMINARy «4 Ch ^ o f F ’ 1 * 4 ' £12.4/ 

SOLUTIONS' ,// 3 (El 3.4 f | 



OF THB 
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SURDQUTfNF RT7(f=SM4| ♦ A L 7P T f NA L l , H I 
COMPLEX FSMAl, A17PT{A2) ,7 
p = 3 . 1 A 1 5° 

TF ( ft F 4 L ( F S^AL ) .GF • 0. ^ ) GO TO 21 
7 = CMPLX( Ar^AG(FSMAL>,-REALfFSMAU ) 
CALL NRAPH ( AL » 7 f Al 2»T t 1 ) . I ) 

N = NAL 7 

TF (H ,FQ, 2,0) N = (NM)/? 

I H .F0« 1 ,0) on T0 JO 

7 = CMPLX( ATMAGC^SMAL » t -PFAL fF$MALM 
can NPAPH f FSMAL »7 ,al/QT(?*N + ;M t-H 
10 DO X J = 1 , K 

Z = rMPix(PMJ-.S) ♦ . 05* ,c^) 

COM NRAPH (FSMAL t Z t ALZffT(J^l| # n 
7 = CMPLX(P*j -, n F*.G 5) 

CALL NPAPH l FSMAl * 7 ,AL7R T< J+N+l) ,11 
^ (H ,F0, l.a) r,o m 3 
7 = CMPLX(P*< J-.5>- .0*, .05) 

CALI NPAPH <FSMAl,7 f AL7RT|J+2*N*2) -i) 
7 * CMPLXI P*J + . P=5 t ,OF| 

CALI NPAPH (FSMAI f 7, ALZPT( j4^*N+2» ,-!l 
^ CONTTNUF 


NAL 7 = ?*N ♦ 1 

J p <H .FO. ?.P| NAL 7 = 4*N + ? 
P F T U P N 


21 IF (F .FO* 1.^) NAf. Z = NAL7*2 
01 7 J = 1 , NAL7 
TF (H . FO. l.nj OH TO 20 
l = CMPLXI P*J- # 0^ , , ) 

CALI NRAPH {FSMA( ,7 ,AL7RTf J + NALZ) I 
?0 7. = CHPlX(P*fj-.5)-.r-5 f .n?> 

7 CALL NRAPH CFSMAL,7 f ALZRTCJ|,l) 
fF (H .FO. ?.P) NAL 7 * NALZ*2 
PF TURN 
FNO 
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SUBROUTINE MRT ( ESMAL ♦ 7«T ,P K t H t NflnnT # 0R N f A M f OHS , 7 S # NS ) 
COMPLEX F<mA! t 7PT( 4?) f l SI 4? ) 

CnMPL EX * I IS 7 f F ,P, r TAN7 f C SFC 7 ,RR f HE S 

»FA|.«4 CRN ( 4 ? 1 * PR ( 4 ? } 

OIMTNS ION N N ( 4 2 ) f M S I 4 .? ) 

RE At *S 7 P , Z T * S ( ? ) . 0 M A P -l*DAM»A»Rn ? TMM # ALtR 
FOU IV At PNC F I 7 « Ml ) I 

n^s = no. 

A = OHS*2.FA/H 

00 = PK 

OP S = F S.MAl 
LH = H 

NM AC. H = A R S ( A M ) / . n S 4- 1 , " 

DAM = AM/NMACH 
NO = 0 

01 IBP l * } , NO nn t 

OR (L» = p.p 

I = 7PTIL I 

s i c = sin 

if ( sic .of. i.F*>p) r,n in i op 

TMM = l.pn 

7 R r ( l I = fmodh i.F?r f \oi 

omath = o.oo 

N = 1 

TF I LM H .or. NPOOT) N = -1 
IP I AM e fO. p. f' I on Tn ior 

00 RO J =1 » NMACM 

DM AC H = CM ACM 4- nAM 

1 = P 

TMM = l.D* - HMACH^OMACH 

on or K=i # i An 
zo = sm 
7! = sm 

IF ( 0 A 8 S { Z° I .GT. n.4010 .PR. DARSCZn * GT . 1 7 4. DP l Gn TO 180 
IF (N .GT. P) C.TANZ = -CDS! N( 71 /COCOS! 71 

IF IN .IT. P) CTAN7 *C0nSI 7 1 /r.OSTNC 7) 

fP (N .GT. PI C SF r 7 = l./C0CPS(7) 

IF fN .IT, P) C 5 F G 7 = !./rnSIN(ZI 

BO = rCSOOTI l .n^-Tv»M*Z*Z/(PD*00*A*A| I 
F = ' OFS/FMW/TMM* n .-OMATH^RR ) **? 4- 7*CTAN7 
p = R0*PC*A*A*T M M 

p * ?.r)p*DFs*( l .n r - pm a ch*rp i*r>MArH/Po*7/o 
P = p ♦ f T AN 7 - 7*C SFC 7 *C SFC 7 

carp = rcARsrp i 

IF (CARP .IT. l.F-40) GO TO IBP 

7 = Z - F/P*.7QP 

IF IS(l) .GF. p.PG .AND. sm . Gf , P.on GO TO } S 
1 = 1 + 1 

IF (I .GT. n r,n TO 1 pp 
z = Z ♦ F /P * .S[>P 

IF IS(I) .LT. P.P^ .no. SI?) .L T . O.OP) GO TP 18 " 
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*1 • tT * UD ‘ 6 • AN0 * tl*BSI2!-SC2M .LT. 1.D-6K- 

GO TO IBP 
q 0 OMTfNUF 
^ « sm/i 

7 I = S( ?)/A 

AL = R 0 *R P - TMMft(7p*7t* . 11 * 71 ) 

NR = NR ♦ 1 

7R TCNR ) *\ ^ PnC ^ nS ° R T( * 5 0P*DABS COSQRT f AL*Ai >B*R j -A| | )/TMM*2,c., 

NN(NR| r M 
190 CONTINUE 

on *5 0 J = l f NP 

on 40 I * l t N» 

lF C OB f I ) ,L T. OR S J I s * r 

40 IF ( OB ( I ) *LT * ORS) PBS = f)B { 1 1 

OBN( J| = DBS 
OBS « 1.B1P 
ZS< J ) * ZRT< TS) 

N^fJ) = NN ( I S ) 

*0 DBMS) * 1 ,F?B 
NRnnt = np 


PF TURN 
FNO 
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c 


SUBROUTINE N RAPH { F SMA t t 7 St ZF *N ) 
COMPLEX*? FSMAL *7S,7F 
COMPLEX* If 7 *F ♦ ° »c T AN 7 t C SEC Z 
RE AL * 5 « 7 » 7 I *H { ?) 

FOUIVAL FNCF C ? t H < 1* ) 


R? 

?! 

IF 


7 = 7 S 
J * 0 

nn * o t * i # i o n 
= HU) 

= H f 2 ) niaci 7 t i FT 174.670^) GO TO 77 

I0ABSIR7i .07. 7.57710+15 .OR* OABMZI) * r »T. 1 

if n iVtInz = - ;psiNin/r.ncns«n 

!F In *.lL 0 1 C TAN7 = ;o;os(7)/cdstn«7) 

IF ( fsi *G7. n ! C5FC7 = l*/CDCO$fZ) 

TF (N .17. CSEC7 = l./COStNIZ) 

F = FSMAL ♦ Z *f 7 AN 7 
p = fTANZ - 7*CSFCZ*CSFCZ 

|F (COABSIP! • L 7 . 1 .17-40) GO 70 72 
1 = 7 - F/P 

!F < H( 1 I .OF. 0.0*" .AMO. H<?l -OF 
J = J + l 

| F <J .07. 6 ) GP 70 2? 

(V,HM, F ^:ToCO .op. HI 2. .LT. . r ; nco. GO TO 72 
b IC|P4BS(P7-HUn.l T. 1.0-6. »N0.0ABS(7!-H<?>». ' T - ’ & 

0 CO NT IN OF 


n.rop) on in 


? ? 7 F « CMPLX ( 1 « F 7" f ? . 0 ) 

PFTiJPN 

If I H( 1) .17. l.O-P .AMO. H ( ? ) 

RF TURN 

FNO 


.IT. 1.0-°) ?F 


CMPtxt i.F 20 ,r- 0 ) 


REPRODUCIBILITY OF THE 
ORIGINAL PAGE IS POOR 
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CHMOLEX 'fsmaUR^? 0 ,nfl, ' V ' 4M * F, ’0«C,M f nHS,DHH,FSMAUH t F<;» 

«f r.f \Tn'<;7; r IT’ c$FCz ' n *’ df s 

FOlirVAlFNCE iJ.sini ,K,P, ’ YK ' T '’"."AM, Al.B 

A = DHS 

P j * 3.141 ->( 1 , 

R < ~ 2 . 00*0 f* pp Q/p 

DAM = AM 

V< = t> l*f M - 1 ) /OnH/^.54 
DFS S FSMAL 
DR = 0.0 
7. = R 

.GF. !.0?0) 3 F TURN 
p - CMP| X( 1 4 f?o f n ./) j 

tmm = i .or - dam*dam 

(am • FQ # n.o) r,o to irr< 

J = 0 

DO 91 k = 1 f i r»rj 
*7 = S(l| 

7 l * %t ?. ) 

!: S :K Bl :Si 

“ COSQRTf !.no_TMM/OK/PK*f 7*j/a /A 

* yK * VK ” 

F = F ♦ 7*CTAN7 

0 = RK*PK*A*A*TMM 

» = p*^7TlN7 n 'r' nAM * BR ’* 0iM/RR * 7/P 
ir irnlvV'ty ~ ? *csec2*ns ec? 

1 .‘«ri :• * lt * ! ‘ o_4 °> rptusm 

(SC1> . GF . r '.Df' .AND cf ?» rr 
J = J ♦ 1 NL * S * 2 * * GE * D.Ori r ,0 TO 16 

f p ( J • GT . 6 1 P FTIP N 

7 * ? ♦ F/P * . Sftr 

fF < SM I .LT, 0.0^ 9 p c/o* , - ^ 

1* fF ( DA R S(P 7- S( 1 ) j IT "i *n S i ? . ;.i T * n * D0) RTTURN 

<?! CONTINUE ' -l-D-6 .AND. DABS( 7! -S (2 )) . | T. 1 . 0 _ 6 ) G n T 

rfturn 

IDO ZP = sc 11/4 

= sr?>/4 

R = z 
FND 
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Si JRRQHT INF POPOEP f A*, a k , G , MC UX 9 7 , NN ,PP ) 
chmolf x*p ? f ?n, 7 K<?n , ;np ( ?n , zknr ( 2 i) ,71 
DIMENSION NN( 1 * ) ,*W( 14) t DP (1 *) 

WRITE ( ) 

DO 1? ! = l, M 0U* 

71 = AK*AK -n ,-AM*AM)* l( I) *Z U 5 /G/G 
\r\ 7K(I) « - A M * AK * f. SO* H Z1 I 

no in j = l y mc u x 
RIG - - 1 *F 2 ° 
nn 70 k - i.mcux 

IF (PEAL(7K(K)> .IT. RIG) GO TO 20 
RT G = PPM ( IV [ K ) ) 

*S * K 
^0 CONTINUE 

7NP( J ) = 7 (KS) 

7K N» (J ) = ZK< X S ) 

DB(M = -7, 54* «.6«*AIMAGI ZK ( KS) ) 

NNB< J ) = N N ( K S ) 

■»0 Z< ( K S ) = C^PLX(- 1 ,F?B,0. 0) 

DO 40 L * l,nru* 

7A L) = ZMP(L) 

/Kill = 7*NP<tl 
NN(l) = hNR(t) 

40 WRITE (f,?I 7(1 ),ZK(L ) , OB (l ) 

PPT!|RN 

1 FORMAT ( * 1 7 , K 7 ,00 S LIB ROUTT NE PC«PFP‘) 

1 FOP MAT ( « n * , 1 p *F 1 4.5 ) 

FND 



SUBPOUT TNF RORDER < AH* A K , G , MC UX « Z f NN 9 OB J 
COMPLEX *8 Z( 21 >,*K{?1 ! V /NR (211 « 2KNR <211 t?l 
Of ME NS TON NN( ! 4 I *NNR ( 14 ) *06(141 
WR |TE { 6,11 
00 10 I = 1 1 MC U X 

ll = AK*AK -(l,-AM*AM|*Zm*ZIII/G/G 
10 ZK < T > ={-AM*AK + C SQR T< Zl ) ) / < 1« 0 - A M*AM ) 

DO ^0 J * l.MCUX 
BIG = -1.E20 
00 20 K = l,MCUX 

TF (REAL(ZK(K) ) .IT. BIG) GO TO 20 
BIG = PEAM/.MK)) 

K S = K 
20 CONTINUE 

ZNR( J) = Z(KS) 

ZKNP(J) = ZKIKS) 

DB ( J I = -2 . *4* 8.68*AIMAG< ZK < KSII 
NNRI J I = NNIKS) 

30 ZK(KS)= CMPLXf- 1 .E2 5 f 0. 0) 

DO AO L * l.MCUX 
Z < L 1 = ZNRCL) 

ZML) * ZKNR(L) 

NN ( L ) = NNR (L) 

AO WRITE (t»3) ZIL ) , ZK (t ) ,DB (L ) 

RETURN 

l FORMAT <M Z.KZ.nB FROM SUBROUTINE P ORDER 1 ) 
3 FORMAT { *C • , 1P5E14.5) 

END 
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TWO WALLS LINEU POLYIMIDE 


99 

8 

3 

2 .C OOOOO 
135*0 OOCOO 


DUMBER OF SPLIS 

MAXIMUM NUMBER OF SOFT WALL MODES CONSIDERED 
MAXIMUM NUMBER OF HARO WALL MOOES CONSIDERED 

NUMBER OF WALLS LINED 

UNIFORM SPL 


l.COCOQD RMS PARTICLE VELOCITY 

3.3 REORDER OPTION 

3.0 OVERRIDE MODE OPTION 


12.CGLJ00 TOTAL DUCT LENGTH 

<►.030000 INCREMENT LENGTH 


12.C3COQO HARD WALL DUCT HEIGHT 


1*400000 GAS CONSTANT 

5 2 3 .C GCGOQ TOTAL TEMPERATURE 

15.C0G0G0 TOTAL PRESSURE 


2003.COOUO0 LOWEST FREQUENCY CONSIDERED 

40C )• CwtOOU HIGHEST FREQUENCY CONSIDERED 


1 *L DO 000 
0.0 
0.3 

6.CCOOOO 


0.530003 


MODEL SELECTOR 
MACH NUMBER 
SOFT WALL OUCT HEIGHT 
NUMBER OF POLYIMID PLYS 
BACKING SPACE 


BEPRODUCIBILITY of the 
ORIGINAL page is poor 
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*4Ch Nu^atas 

0*190000 0*200000 

SOFT MLL CUCT HEIGHTS 
5*799999 6.000000 


0.220000 

6.200000 
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DB DIFFERENCE AFTER 4.00 INCHES OF LINING 


2000. 

2.60 

2100. 

4.09 

2200. 

4.30 

2300. 

4.60 

2400. 

4.98 

2500. 

5.43 

2 600. 

5.92 

2700. 

6.46 

2800. 

7.04 

2900. 

7.63 

3000. 

8.15 

3100. 

8.91 

32CO. 

8.53 

3300. 

8.21 

3400. 

7.93 

3500. 

7.65 

3 6 3 J • 

7.33 

37C0. 

6.58 

3800. 

6.62 

3903. 

6.28 

4000. 

5.57 


□ B DIFFERENCE AFTER 8.00 INCHES OF LINING 


2000. 

4.93 

2100. 

7.10 

2200. 

7.78 

2300. 

8.55 

2400. 

9.44 

2500. 

10.41 

2600. 

11.46 

2 700. 

12.54 

2800., 

13.64 

2 9C 3 * 

14.72 

3000. 

14.88 

3100. 

14. C2 

3200. 

12.60 

3300. 

11.59 

34C 0. 

10.82 

3500. 

10. 21 

3600. 

9.68 

3700. 

9.20 

3800. 

8.76 

3900. 

8.36 

4tiOQ . 

7.99 
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UB 0 1 FFeRfcNCE AFTER 12*00 INCHES OF LINING 


z 000 . 

7.03 

2100. 

9.68 

2200. 

10.80 

2300. 

12.05 

2400. 

13.43 

2500. 

14.94 

2600. 

16.52 

2700. 

18.04 

2800. 

19.53 

2900. 

20.79 

3000. 

19.89 

3100. 

17.62 

3200. 

15.48 

3300. 

14.01 

3400. 

12.91 

3500. 

12.05 

3600. 

11.35 

3700. 

10.76 

3800. 

10.24 

3900. 

9.78 

4000. 

9.36 


FREQUEN 

CIES 

FINAL SPECTRUM 



2000 

ro 

2300 

127.97 

125. 32 

124.20 

1 22 ;95 

2400 

ru 

2700 

121.57 

120.06 

118.48 

116.96 

2800 

to 

3100 

115.47 

114.21 

115.11 

117.38 

3200 

TO 

3500 

119.52 

120.99 

122.09 

122.95 

3 600 

TO 

3900 

123.65 

124.24 

124.76 

125.22 

4000 

TQ 

4000 

125.64 




OB DIFFERENCES 

I SMOOTHED) 

AFTER 12.00 

INCHES OF 

LINING 


a **5 10.24 11.42 

15.73 17.28 18.79 

18. 75 16.55 14. 75 

11*70 11.05 10.50 


12.74 14.18 
20.16 20.34 
13.46 12.48 
10.01 9.57 


m 



APPENDIX IV 


PROGRAM FLOWCHART 
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LL 







© 



reproducibility of the 

ORIGINAL PAGE IS POOR 
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7 


& 


\ 


fliiRmjut jhpqo ir. t, 
iric , i spl. in 

rn — 


C f WEWntC - PWWRItll SHE ET 


Will Tp 


CIi IXWJJ nti^sjH rxi 
: P*3,2E-4«10.»« tfiBS I17L1/3IM 
W5-S.^]-F/C-ISr-2,SM 

n*»s 

IMD 



tn.tt/sis. 

rp 2 *Pi/ 2 ./JMn a.i 

ItL^i/w.T 

E» 1 .fl25l« U./PBR— 2 -tl VtXP l.50T2»Trtl/0Un 
AVOC°.OT1«'!m,'BEL/’PM""Irt«rTM/ iTm*.416) 

RV2^flVDC/2, 

ilSlSo • ,v? ' 3an ' 1 »wa«« 2 H. 125 *- , BSS'tJFV 

Z • D .5 



, SN « 2.SW « Olfl « F « PflA/VP 
| I1AC*PP2«Vp/C«E.EXf 

' n£-5»Trt.«.lS«S8flT IFU nwvOJBM.-PBHW PBfl.JJftl IJQ.« |1H 

I irz»mwr»fWAC»flVDCi*<iBr 
n.T».< 5 «B]B^ n.-.i-nm p»u-». 2 C-s/pm«vp 



T 

IB TKfl*2.«PJ«F/[«B3P*2,54 

rz«CAPu< mi/BOC.nx/BBC-CDT nxot > 
return 

FI1PV - BJRC/Vp«U.*T.B«3»«5IO«fla: ! 

wpv»-w n**m-n«».jE-5/FJfl i 


l 



ZI-CKPIX IPI,PX-BIC«CBT n«CI l 

me » cabs i?n 

SBPTY - SCAT IVPCFwJ f P/ Jflji «« 2 ) 
FWSoVP- SJBTV 
fPBI * - 2 .i'BA&«*> 4 « KT«BTPVtfU«RXPVI i 
F*P - 1,3 - PP8I..5/5SB1V 
»pK-yp-rip/Fi«p 





III. tt/RSC 

rr»vp» 

RETURN 





< wne i». in j». jin > 


n - ii/Btc 

RETURN 
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\_g j _ gflauu jg UX IF. flu . M. c, Dob. nCux. ‘Vj. *j. j, fMt, n — 7 

— X ~ 1 


L f «««« an n Kiti/t 




si » cnPLX u.fl.i.n 


: W-2.-Pi.F/'C 
! 5HS * S 
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ajafiairriwc i-wobs w.oi.ii.m.fNi.c.h. h fji 


t R • OMW.SM/h 
: fi • i. iKisjjssit 
: g • 2.UM»j.Ffii/r 

1 If** p ^’ W-l>/IMti/2.SK 
ffl « FShAL 




5 - fNPtKii.no. o.oi 

1MM w v.OO - OfiN«Df»i 


A1 • Sill 

II • 0 121 


<C . ra SI2. # !H . - C1 - 3-5111015 ,0 

\5 - mis ran .6i. m. nooi 


[ reiIkin | 



hr*? - -cbsjk m /cocos m 


n«K2 .cocos m/cosjN m 


I cscci - i, /cocos m 


CSEC1 • l./COSJR 111 


J** ij!5?5 T n . rs-ihn/hK/Ritw rjwi/p/Fi ♦ ik.ykii 

s-'tnn/ihnw n.oo - awwBfliwwj 

f * f ♦ 2 «nntr? 
t • FKwWlwfWftwTtW 

J * 2. OOwOF J« >1. DO- Wh«8R) wDArt/lfWI/P 


| <3(M»5IP1 .LI. 1.0-KJJ. 

1 

: p_ 


S 
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